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ABSTRACT 


This study develops simulation models for the evaluation of the traffic performance 
at a three way intersection with no signal. In particular, the models were designed to test 
and to evaluate the wartime road capacity of the current and a proposed road system in 
the Korean theater. The models describe how traffic at the road intersection will vary 
with changes to the input variables and priority rules. An analysis is performed to com- 
pare performance measures of the intersection given various traffic conditions. Specif- 
ically, the analysis will enable an efficient modification of the roads or an appropriate 


road usage plan in wartime. 
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I. INTRODUCTION 


A. BACKGROUND 

The wartime traffic on the roads in the southern theater of the Republic of Korea 
needs special attention, especially the traffic involved in a war which is very dynamic, 
and in particular a war which involves a sudden attack from the north and the necessity 
of falling back to the south. There are several main avenues leading from the front, from 
which many streets to the east or west are derived, so we have many three way inter- 
sections with undeveloped roads. The traffic on these roads is not very heavy in peace 
time, but it appears as if the limited road capacity may restrict the trmely movements 
of forces in a time of emergency. 

The study of traffic in this specific area has never been done before because under 
the current strategy no forces are allowed to withdraw, with the exception of supply 
forces and forces needing to be replaced. We assume that in many cases all the forces 
will have to move simultaneously or in a certain sequence. For instance, the enemy may 
select battlefields which are quite different from those we expect. Also, the deployed 
forces may move forward or backward simultaneously as the situation changes. It is 
strongly believed that traffic congestion will occur in the assumed contingencies, and 
therefore, it 1s necessary to have some idea of the dynamics involved. A particular 
question which could be asked is, under what condition will the traffic moving onto the 
main road from the side road be unstable. That is, when will the queue of traffic continue 
to grow in an unstable manner. This thesis presents a very limited approach to the 
problem of the dynamics of force withdrawal from the front or a deployment of reserve 
forces to the front, but it 1s an important first step to the understanding and solution 


of the whole problem of traffic control and flow in the Korean theater of operations. 


B. DESCRIPTION OF THE PROBLEM 

The purpose of building a stochastic simulation model is to describe traffic per- 
formance at a three way intersection. This type of intersection is a principal capacity 
limiter and main source of delay of the current road system in the Korean theater. At 
first, an analytical model for the current road system under special simplifying assump- 


tions 1s developed, including a transition diagram and brief solution procedure. We will 


assume that a unit, which possibly consists of many vehicles, 1s treated as a single cus- 
tomer as in general queueing problems. Each unit’s transition time is equated with ser- 
vice time in the sense of a classical queueing situation, |Refs. 1,2: pp. 13-20|) 9 ie 
interarrival and service times of a unit are assumed to have certain respective distrib- 
utions, which will be discussed in detail later. Further, the interaction of different units 
at the intersection will be controlled by traffic policies to be described in Chapter II. 

Due to difficulties encountered in the analytical approach, two simulation models 
are developed for the current and a proposed road system. We are concerned with 
finding how the traffic performance at the intersection varies with changes to the arrival 
and service time distributions. We consider performance in terms of the average waiting 
time of the units and the average number of units waiting 1n each direction of the inter- 
section. 

The intersection system must be thoroughly understood because blockages at an 
intersection may cause serious problems in conducting operations in an area where units 
are deploved with high density relative to the road capacity. The units moving or wait- 
ing because of blockages are very vulnurable to enenyy ambushes since they are not fully 


ready for engagements. 


C. SCOPE OF THE THESIS 

In attempting to achieve the goals of this study, an analytical model is developed 
up to the transition diagram, together with the general solution procedure outlined in 
Chapter II. It is found that the simulation model is needed to solve the problems be- 
cause of the difficulties in handling various inputs and because of the limited solution 
capability available for the analytical model, 

Two simulation models are built with SEMSCRIPT IJ.5 in order to conduct the 
simulation experiments. The first one is based on the current road system and the other 
is for a proposed modification to the system. The two models are presented in Chapter 
Iif and Chapter [V. Figure 1 and Figure 15 depict the current and proposed road sys- 
tems respectively, with possible movement routes. 

Chapter V will analyze and compare the results of the two simulation models. The 
analysis is done by exporting output files of SIMSCRIPT runs to the GRAFSTAT 
program on the IBM 3033AP main frame. This allows one to graphically examine the 


models for stability of the performance measurements and to examine the results in the 


stable cases as a function of several input factors. Finally, conclusions and recommen- 
dations are presented in Chapter VI. 


II. PROBLEM STATEMENT 


As mentioned before, we can apply this model to the wartime traffic in a theater of 
Operations where the road capacity is limited 1n comparison with the potential road us- 
ers. We assume that there are no traffic signals and no military police for traffic control 
and assume services are controlled by service policies we set up for the simulation study. 
It is the fundamental road capacity and traffic performance at the intersection that we 
want to examine, therefore, we must fully understand the following assumptions to 
properly grasp the specific features and logic of the model. A current intersection 1s 


presented in Figure 1. 


A. ASSUMPTIONS 
1. Units 
A unit consists of a group of vehicles, the number of which depends on the or- 
ganization of forces and the tactics of a war. All the vehicles in the unit will be treated 
as one customer at the intersection (general queueing facility) because the movement of 
a combat unit such as a company or a battery must be conducted under tight command 
and control to minimize confusion and loss of vehicles. Uhus we assume that no unit 
or vehicle can intervene between another unit’s vehicles. 
2. Arrivals and Queues 
Poisson arrivals are assumed for the analytical model. Each arrival from the 
three directions is filed in an appropriate queue if it is not its turn for service. Our 
intersection queueing facility has three queues; a north, south and east queue and there 
is no interaction between vehicles in each of the queues, that 1s, a First Come First 
Served principle is effective within each direction. Actually the time of arrival is that of 
the unit to the end of the line, potentially some distance away from the intersection. 
Since units move up to the intersection in close order as units which arrived previously 
are served, the distance and stop and start effects of the queue at the intersection are 
ignored. In fact, the waiting time of a unit is proportional to its distance from the 
intersection. This could be exploited in the analysis but will be ignored. 
3. Route Selection | 
Possible routes for a unit from each direction are shown by arrows in Figure 1. 
One of two routes will be selected for a unit from each direction and this selection is 


made when a unit becomes the first unit in the queue and no unit is presently served 


North route] 


East route2 
East routel 
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North route2 
South route2 


South route 


Figure 1. A current Intersection.: Traffic coming from the north either contin- 
ues south or turns east if it is not blocked. Similarly traffic from the 
south can continue north or turn east if it is not blocked and traffic from 
the east turns either north or south. Priorities for service on the routes 


will be discussed in detail in Chapter II. 


from the direction or when an arriving unit sees no unit in the queue. The probability 
used to decide a route 1s an input variable in the model and a uniform (0,1) random 
variable L will be drawn to determine the route selection. If U is less than the proba- 


bilitv P of taking route], route! is selected, otherwise route2 is selected. The probabilities 


of taking north routel, south routel, and east route! may be different from each other 
and are uSer input paranieters. 
4. Service 

Service begins when the first vehicle of a unit starts from the intersection and 
ends upon the passage of the last vehicle. Once served a unit is assumed to leave the 
system, even though, in reality, it may significantly affect the service conditions of follow 
on units. Service umes depend on many factors, such as the length of a unit, road con- 
ditions, weather, number of refugees and so on. Though many of the factors should be 
studied separatelv, they are assumed to be compounded into a single random variable 
used for service time in this model. Our main interest in this thesis is the interactions 
between units at the intersection. We assume the service times for the six routes are 
independent and exponentially distributed in the analytical model. This is done for an- 
alvtical tractabilitv and is probably not physicallv valid. More general and probably 
more realistic assumptions are made in the simulation models. 

5. Service Policies for Model 1 

In Model I, service for a unit starts immediately upon arrival when there is no 
blocking to the selected route, otherwise it 1s filed in the queue. There are two possible 
routes from each direction named route! and route2 as shown by arrows in Figure I. 
Suppose that the set (north route, south route, east route) indicates the state of move- 
ments at the intersection where each entity within that set can have value QO, | or 2 ac- 
cording to the routes being used. Here 0 indicates that neither of the two possible routes 
from each direction are used; otherwise | and 2 represent the route number used. It is 
clear that all combinations of 3° (= 27) states are not included in the possible movement 
states of the model because of the interactions between the routes. For instance, the 
state (1,1,2) never happens since south routel and east route2 can not be used simul- 
taneously. 

In Model I, 1n fact there are onlv 14 states of movement which are acceptable: 
(0,0,0), (0,0,1), (0,0,2), (0,1,0), (0,2,0), (0,2,1), (0,2,2), (1,0,0), (1,0,2), (1,1;0)) (aes 
(1,2,2), (2,0,0), and (2,0,2). Even though these sets are not themselves enough to study 
the quantities of interest, the service polices are based on these possible movement 


combinations. 


A movement on a certain route may block more than one route to which units 
have to move, therefore the priorities for movement must be specified according to the 
route for which a service is just completed. Note that we have assumed there is no traffic 
signals and no military police for traffic control. A unit desiring to move in the system 
must walt until its route 1s clear and it has priority (in the case of ties). We break ties 
in cases in Which more than one unit has conflicting priorities simultaneously. The 
routes that block movement and the next priority for a service given a blocking unit 


which has completed service are shown in Table 1. 


Table 1. PRIORITIES AT THE END OF SERVICE OF A BLOCKING ROUTE 
(MODEL 1) 


Route which 1s about on ; 
Next Priority, in order, at the end of blocking 
to finish service 


north routel east routel, (north routel,2) 
north route2 east routel, (south routel,2), (north routel,2) 


south routel (east routel,2), north route2, (south routel,2) 


south route2 north route2, (south routel,2) 


east routel (north routel,2), south routel, (east routel,2) 


east route2 south routel, (east routel,2) 





The left column. of each row in Table 1 indicates the routes on which service 
was completed and the right column of each row identifies all possible routes which can 
be selected by units Waiting to pass when service was completed by the route denoted in 
the left column. Priorities for next service are given to the right column routes in the 
order listed. We set up the specific priorities for this model by breaking ties and placing 
the routes for which service was just completed last in order of priority. Next priorities 
at the end of service on each route will be discussed in detail below. The routes in the 
parenthesis denote the two possible candidates from the direction blocked, so only one 
of the two routes can be Selected for the first unit in the queue. Table 1 also shows that 


a newly arriving unit desiring to move on a route of a left column and row can only 


move if all the right column routes of the same row are not used, otherwise it has to wait 
for service until next, priority 15 given tovit. 

After service for a unit on a route in the left column of Table | is completed, the 
priority for next service 1s conditioned on other movements at the intersection. As an 
example, the third row could represent units that have to move to (east routel or east 
route2), north route2, and (south routel or south route2) and are blocked by a unit 
presently using the left column south routel. Upon service completion on south routel 
the first check 1s to determine if there is blocking of east routel or east route2. Let’s 
assume that a unit from the east has been waiting for movement to the north (east 
route2 in the first position of the right column), the unit can move immediately because 
east route2 is always clear right after the service completion of south routel. Next, 
service for north route2 (second position in order) will be considered and a unit waiting 
to go from the north to the east (north route2) can also move immediately because there 
is no blocking to that route at the end of south routel’s service regardless of the service 
on east route2. Note that east route2, which has the first priority for service, does not 
block north route2. Finally, service on the south routes (third position in priority order) 
will be considered and movement to either south routel or south route2 will be condi- 
tioned on the status of east route2 and north route2. If we had assumed east route! (the 
other route in the first parenthesis of the right column in the third row) was waiting to 
move instead of east route2, we must check to determine if there is a unit on north 
route! upon the service completion of south routel. If north routel is active this route 
also blocks east routel. If there is no unit on north routel, east routel can be used, 
otherwise service on the next route in the order of priority, north route2, must be con- 
sidered. In other words, even though a unit has priority upon a route’s service com- 
pletion, the unit loses priority if it is blocked by route presently in service. 

In Table 2 through Table 7 we tabulate cases of the priority problems. The 
triple in the row designator of each Table describes the routes selected by units which 
are specifically blocked by the route that 1s about to finish service. The first element 
refers to the traffic from the north, the second refers to traffic from the south and the 
third refers to traffic from the east. In each case the element can have value 0, | or 2, 
where 0 indicates no traffic waiting from that direction, | indicates traffic waiting for 
route! and 2 indicates that traffic is waiting for route2. The triples in the column des- 
ignator represent all possible movement states when the route of interest, which 1s 


specified in the Table title, is about to finish service. Finally, N1, N2, $1, $2, El, and 


E2 represent north routel,2, south route 1,2, and east route 1,2 respectively. A route 
identified at the intersection of a row and column has priority upon the service com- 
pletion of the specified unit in the title of the Table. N/A indicates that the combination 
of the row and column condition is not feasible in this Table. A 0 at the intersection of 
a row and column indicates that no additional route can be used immediately after ser- 


vice completion of the unit specified in the Table title. 


Table 2. SERVICE PRIORITIES AT THE END OF A NORTH UNIT MOVING 
SOUTH (denoted by north routel) 





In the last row of Table 2 the entry (2 0 1) 1n the left colunin means that when service 
on north routel 1s completed, north route2 1s selected for the next north unit (the first 
2), there is no one waiting from the south direction (the second 0) and someone is 
Walting to go from the east to the south (the last 1). The priorities for movement of 
Waiting units must be conditioned on the present situation at the intersection, and 
this is the column indicator. Thus while a unit from the north 1s moving south, a unit 
from the south can be moving north (south route], designated by the second | in the 
header of column 3) if no unit from the east is turning north (east route2), and can 
move east at anv time (south route2). Therefore there are only five movement states 
possible when north routel is about to finish service. The priorities for next service 
given the row entrv are dependent on which one of the five movement states describes 
the intersection Just prior to the completion of north routel. Not all row and column 
combinations are feasible and are designated by N/A. Thus in row (2 0 1) and column 
(1 0 2), east route2 is being used so no route 1s selected for an east unit (the 1 of (2 
0 1) can't be selected). A O indicates that no additional route can be used after service 
completion of north routel. Thus for row (2 00) and column (1 2 2), the unit from 
the north can not move east because a unit from the south ts moving east (2 in the 


second position of the column heading). 
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Table 3. SERVICE PRIORITIES AT THE END OF A NORTH UNIT TURNING 
EAST (denoted by north route2): Explanation of Table arrangement is 
provided on pages 8 and 9. 


While a unit from the north is turning 

east {north route2), only east route2 is 

not blocked by the unit. Thus we have 

S1 only two possible movement combina- 

iss tions at the intersection just prior to 
S? north route2 finishing service. All com- 
ean | N/A binations of the 6 routes except east 
route2 are in the row designators, which 
indicate the possible routes selected bv 
E units waiting to pass when service on 
Nt and S!] NI north route2 is completed. A unit on 
By N/A east route2 (the last 2 in column (2 0 2)) 

N]1 and S2 N1 and $2 | still blocks a unit waiting for service on 
Enea > N'A south route! when north route2’s service 
X 2 is completed, thus south routel loses 

4 priority, which is designated by V. The 

e — case of VV in row (2 | 0) and column (2 

S1 N2 Q 2) indicates that the priority initially 

El NIA given to south routel is transferred to 
a2 S2 north route2 because of the unit on east 


El and S82 N/A route2. 


NI Nl 
El N/A 
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Table 4. SERVICE PRIORITIES AT THE END OF A SOUTH UNIT MOVING 
NORTH (denoted by south routel): Explanation of Table arrangement 


is provided on pages & and 9. 
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Table 5. SERVICE PRIORITIES AT THE END OF A SOUTH UNIT TURNING 
EAST (denoted by south route?) 


Explanation of Table arrangement is provided on pages 8 and 9. A unit on south 


route2 blocks only a unit from north turning east (north route2) and units from the 


south routes themselves. Therefore north route2 has the priority for next service if 


there 1s no unit moving on east routel, else the priority will again be transferred to 
the south routes. This 1s designated by V in row (2 2 0) and column (0 2 1). Note 
that in Table 1, north route2 is placed ahead of the south routes in order of priority 
at the end of service on south route2. This is because we have assumed that once 
service 1s completed in one direction, priority for next service will be given to the 
other directions which have been blocked by the unit leaving the system. A unit on 
the east routes, the last 1 or 2 of the triples in columns 2, 3, and 5 allows no service 
On south routel and the movement on east routel (1 of column (0 2 1)) at the end 
of south route2’s service also blocks north route2 (see VV). The 0’s in rows 1, 3, and 
4 and columns 2, 3, and 5 are because units on east routel,2 continue to block service 


of units attempting to move at the completion of south route2’s service. 





Table 6. SERVICE PRIORITIES AT THE END OF AN EAST UNIT TURNING 
SOUTH (denoted by east routel); Explanation of Table arrangement is 
provided on pages 8 and 9. 
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Table 7. SERVICE PRIORITIES AT THE END OF AN EAST UNIT TURNING 
NORTH (denoted by east route2) 


Explanation of Table arrangement 1s provided on pages 8 and 9. A unit in service 


On east route2 blocks only south routel and east routel,2. Therefore south routel 
has the priority for next service if there is no unit in service on north route2, else the 
priority will be again transferred to the east routes. Note that in Table 1, south route] 
places ahead of the east routes in order of priority at the end of service on south 
route2 because we have assumed that prioritv for next service is given to the other 
directions once service from one direction is completed. The unit on north route2, 
the first 2 of column (2 0 2) blocks a unit waiting for service on east route! and south 


routel, which 1s designated by V. 





B. ANALYTICAL MODEL 
1. System Notations 
In this section, we introduce the notations required to develop an analytical 


model. Note that subscripts of ”, s, and e represent north, south and east respectively. 


north units Units from the north direction. 

d, Poisson arrival rate at the intersection from the north direction. 
HL, Exponential service rate of a unit from the north direction. 

O (t) Number of north units waiting for service, as a function of time ¢. 


Next.north.route Next route selected for the first unit in north.queue. 

north.queue Number of units waiting for service from the north and next route 
for the first unit in the north queue. 

north.run Route that a north unit uses now. 

Pe. Probability that a north unit will take north routel. 

Note that north.queue represents two quantities, Q,(t) and next.north.route with 
“!” between them. As an example, “3/1” of north.queue indicates that there are 3 units 
in the queue and the first unit in the queue 1s supposed to use north routel. 

There are also south units, /,, u,, south.queue, south.run, next.south.route, Q, 
and P, for the south direction and east units, /,, u,, east.queue, east.run, next.east.route, 
O, and P, for the east direction. 

2. A Transition Diagram 

Suppose that potential units arrive in accordance with Poisson rate /,, /,, 4, 
from 3 directions and that service rates are exponential y,, u,, and u,. We can then the- 
oretically compute the quantities of interest for the system; the mean number of units 
in the queues and the average amount of time that an arriving unit waits in the queue. 
To accomplish this we must first decide upon an appropriate state space. It 1s clear that 
the state of the system must include more information than merely the number of units 
in each queue of the system and the present movement state (1 of 14 movement states 
we discussed earlier). For instance, to decide the movements after a service completion, 
it would not be enough to specify that there are n,, n,,n, units in each queue and which 
routes are being used, but we would also have to know the routes selected for the first 
unit in the queues. 

To account for this, we define a state variable (north.queue, south.queue, 
east.queue, north.run, south.run, east.run), where the first three entities give information 


about the number of units in the queues and the desired routes for the first units in the 
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queues and the last 3 entities, north.run, south.run and east.run represent the routes that 
are presently being used by north, south, or east units. For example, the following 


sample states have the following interpretation. 


(0/0, 2/1, 1/0, 0, 0, 2): 
The entry 0/0 means that there 1s no unit in the north queue. 


The entry 2/1 means that there are two units in the south queue and south routel is 
selected for the first unit of south queue. 


The entry 1/0 means that there 1s one unit in the east queue. A route 1s not decided 
for the unit because an east unit 1s still being served. 


The last three entries 0, 0, 2 means that only east route2 out of the six routes 1s pres- 
ently being used. Note it blocks south routel. 


(0/0, 2/1, 1/0, 1, 0, 2): 


The conditions are the same as the above state except a new north unit uses north 
routel (the 0 in position four has changed to a one). 


(0/0, 1/0, 1/2, 1, 1, 0) : 


The entry 1/0 means that there is one unit in the south queue. A route is not decided 
for the unit because a south unit 1s still being served. 


The entry 1/2 means that there is a unit in the east queue and east route2 1s selected 
for the unit. 


The entry | in position four means that a north unit uses north routel. 


The entry | in position five means that a south unit uses south routel. 


It should be noted that the process will go to the second state (0/0, 2/1, 1/0, 1, 
0, 2) from the first state (0/0, 2/1, 1/0, 0, 0, 2) when a unit arrives at the intersection from 
the north and takes north route! (which occurs at rate P,/,). Also the process will go 
to the third state (0/0, 1/0, 1/2, 1, 1, 0) from the second state (0/0, 2/1, 1/0, 1, 0, 2) when 
the east unit completes its service and east route2 is selected for the following unit 
(which occurs at rate (1 — P.)y,). 

As discussed earlier, there are 14 combinations of routes which are allowed in 
the system. Anyone of these combinations may have many states since we have ap- 
proximately {3(” + 1)}° States for each combination of the routes where n is the length 
of the queues. We can begin to analyze the states from the following development of 
the (0,0,1) route combination. 
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(0/0, 0/0, 0,0, 0, 0, 1) 
(0/1, 0/1, 0/1, 0, 0, 1) 
(0/2, 0/20/250N@maID 


(1:0, 1/0, 1/0, 0, 0, 1) 
(1:1, 1/1, 1/1, 0, 0, 1) 
(1/2, 1/2, 1/2, 0, 0, 1) 


(n/0, n/0, n/0, 0, 0, 1) 
(n’l, n Pen eee) 
(ni2, 2 Zam 270s Oe 


Of course in this model all {3(” + De states are not possible in the state space 
of the model. For instance, the state (1/1, 1/1, 1/1, 0, 0, 1) 1s not in the state Spacetsimiee 
an east unit 1s still being served, and the direction of the first unit in the east queue has 
not been designated 1.e., the third entry must be 1/0 if a unit from the east is being 
served. Asa result, each route combination may have less than {3(n + 1)}° states but 
this 1s still too many states to describe in a transition diagram for a whole model. In 
addition, it 1s not straightforward to build a total state space since we would have to 
consider all interactions of the routes for each movement state. 

Rather than a direct attempt to solve the whole problem we can observe how 
the transitions occur and what the balance equations look like by starting from the tin1- 
tial point, that 1s, the base state (0/0, 0/0, 0/0, 0, 0, 0). A brief transition diagram that 
describes a small part of the whole model 1s shown in Figure 2. 

3. Solution Procedures 

Up to now we have decided upon the state space and tried to make a transition 
diagram, although they are not fully described. Also we have assumed Poisson arrivals 
and an exponential service rate for the system, so the model can be analyzed as a 
continuous-time Markov chain, [Ref 3: pp. 210-273]. 
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Figure 2. <A partial transition diagram of the model: The short arrows at the 
bottom of a state represent manv more arrows pointing in and and out 
of the state. The notations of Pn, 4m, and un are used for p, . 4,, and 


u,. The notations of Ps, 4s. and us are used for the three parameters of 


the south and Pe. /e. and ye are used for the three parameters of the east 


direction. 


Consider limiting probabilities p, for all states. p,= lim p,{t) where j 1s inde- 
pendent of state { and p, is the probability that the process enters state j from state /. 
If we assume that the limiting probabilities p, exist, then they can be obtained by the 


following formula, [Ref. 4: pp. 246-248]. 
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v; py= »% Pr Py for all state j, k. (2a) 


kj 
Pi, for all statey, where (2 
j 
Vi Py ; Rate at which the process leaves state /. 


>" Pry Py» Rate at which the process enters state /. 
#] 


V : Leaving rate from state /. 

Vp ; Leaving rate from state k. 

Pj : Proportion of time the process is in state /. 

Pr ; Proportion of time the process asamestaren. 

Py ; Probability of transitions from state & to state /. 


Once We solve the balance equations, it is easy to answer our questions in terms 
of these limiting probabilities. Liowever, as we see in Figure 2, the number of states de- 
rived from the base state increases geometrically. Suppose the maximum number of “n” 
is 9, then we have approximately (10 x 3) x 14 = 415,800 states and the same number 
of balance equations where 14 indicates the number of all possible combinations of 
movements. It is not straightforward to develop the transition diagram and to write all 
the balance equations and as a result, the approach to solving the problem necessitates 


the simulation method. 
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It. SIMULATION MODEL 1 


A. ASSUMPTIONS AND NOTATIONS 

The assumptions and notations of simulation Model 1 remain the same as those for 
the analytical model discussed in Chapter II, except that the interarrival and service 
times may now have a Gamma distribution with various parameters, [Ref 5: pp. 
267-273]. 


B. INPUTS TO THE MODEL 

Our primary goal is to examine the waiting time and the number of units in the 
queue of the steady state under certain conditions given by specified input parameters. 
The simulation model is made ready to run by entering the necessary input values. This 
part of the model is extremely important since it provides the user an opportunity to 
specify the particular conditions which are to be examined, as well as the values which 
establish the boundary conditions of the simulation. These input values represent fea- 
tures of the system such as the probability of taking routel and the shape of the service 
and interarrival time distributions. Hence, for Model 1, we have a total of three factors 
which we divide into different levels in Table 8 for the purpose of controlling the size of 
the simulation. 

We assume that successive service and interarrival time distributions for the three 
directions are independent and identically distributed. In addition, we assume the 
probability of taking routel is the same for each direction. Thus from Table 8 we have 
2 levels of @,,k, and &,, 3 levels of P, and 3 levels of traffic intensity. Thus we initially 
simulated (levels of 6,) x (levels of k,) x (levels of k,) x (levels of P,) x (levels of p) = 
2X2x2x3x3 = 72 cases of system conditions. Note that we have 4 different com- 
binations of the Gamma shape parameters k, and k,; (1,1), (1,4), (4,1), (4,4). The value 
1 for k, or k, represents the case where service time or interarrival times are skewed, i.e. 
the arrivals form a Poisson process and the service rate is exponential, as in the analyt- 
ical model. The value 4 for the shape parameter indicates that the Gamma distribution 
is regular, [Ref 6: p. VIII.7]. 
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Table 8. MAJOR FACTORS AND LEVELS OF SIMULATION. 


Service time distributions 
Gamma(@,,k,) where 


(denoted by the subscript | 
Cie 2 55.0 anda es lee 


On parameters @ and k) 


Traffic intensities p 
mean service time ) 


0 = 0.5, 0.4, 0.33 


(p mean interarrival time 


Interarrival time distributions 
Gamma(6,,k,) where 


denoted bv the subscript 2 
: = 0,=— and k,= 1,4 


on parameters @ and hk) 


Probabilities of taking route! P, = 0.25, 0.5, 0.75 





C. MODEL IMPLEMENTATION 
1. Simulation Language 
The model is written in Simscript II.5, which has become the standard simu- 
lation language for military applications and also a powerful and flexible simulation tool 
in all areas of military modeling. Simscript I[.5 reduces simulation programming time 
and cost several fold compared to FORTRAN, and is enormously more efficient and 
flexible than interpretive simulation techniques, [Ref. 7,8]. 
2. Model Structure 
a. General 
The program contains a preamble, a main, 4 routines and 9 event routines. 
The 4 routines include read.data, initialize, single.run.report, and sum.report. The 9 
event routines are north(south,east).unit.generator, new.north(south,east).unit, and 
north(south,east).service. Routines are called by the main routine or by other routines 
and return control to the calling routine. Event routines are special types of subroutines 
which are accessed by being scheduled, {Ref. 9,10]. The functions of each routine and 


interactions among routines are as follows: 


b. Preamble 
It contains the definitions of entities, attributes and sets, and other global 
information. 
c. Nain 
This is the first routine to receive control during the execution of a simu- 
lation run. Its function is to initialize the svstem and to combine the input parameters 
for performing a series of simulations. It also computes and prints some simulation sta- 
tistics. 
d. Read.data 
This routine reads all the input parameters provided by a user. A user 
specifies these parameters in the main routine. The input parameters are three seed 
stream numbers, the desired number of units to be simulated, the desired number of 
replications, the various Gamma distribution parameters for the interarrival and service 
times, and the probability of taking routel for each direction of the intersection. In this 
thesis we examine only the cases described in Table 8. 
e. Initialize 
The spaces for all array elements are allocated. It also schedules initial 
events. 
J. Single.run.report 
All information of a single simulation run is collected here and a verification 
of the model can be dore by examining this report. A sample report is shown in Table 
11 of Appendix A. 
g. Sum.report 
This routine summarizes the output from single.run.report. The report can 
be used to note the stationary statistics and to observe how the output statistics change 
with respect to the input parameters. In addition it can be used to analyze the boundary 
conditions for stability; this being one of the research objectives. The study of this re- 
port is used to determine an appropriate design size for Model 1. A sample section of 
the resulting report is shown in Table 12 of Appendix A. 
h. North(south,east).unit.generator 
Note that south and east in the parenthesis indicate that the same routines 
also exist for the south and east directions. New.north.unit event routine is scheduled 
immediately after the north.unit.generator event routine is called and this event itself is 


rescheduled for a next arrival. The interarrival time for the next north arrival is a 
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Gamma random variate whose parameters, as described previously, are inputs to the 
model. The termination of simulation 1s also controlled by this event routine after the 
desired number of units have been simulated from each of the three directions. 
1. New.north(south,east).unit 
A north unit is created immediately after this event routine is called and this 
unit is either served directly or filed in the queue. Whether it is served or filed in the 
queue and which route is selected for the unit, depends on other movements at the 
intersection and the route selection procedures we discussed in Chapter II (see route 
selection procedure on page 4). If the unit is served immediately, a north.service event 
routine is scheduled to occur at the end of service. Note that the service time is a 
Gamma random variate. 
J.  North(south,east).service 
A unit that has passed the intersection is destroved immediately and this 
unit leaves the system. The units blocked by the leaving unit are served according to the 
Service priorities discussed in Chapter IJ and the appropriate north(south,east).service 
event routines will be scheduled based on new service times (Gamma random variates). 
3. Logic of the Model 

The flow charts of 4 major event routines for the north and a macro flow chart 
for the whole madel are depicted in Figure 3 through Figure 7 respectively. The flow 
charts for other directions are omitted to avoid redundancy. The logic of the event rou- 
tines for the south and east is exactly the same as that for the north, though the condi- 
tions for movement to a given route is quite different. 

We can divide the north(south,east).service event routine into two phases; a 
route selection phase for the first unit in the queue and a decision phase for the next 
movements. There is no difference between the north, south or east directions in the 
route selection phase. However, the decision phase is quite different for each direction 
since each route blocks different routes and the conditions of movement to a blocked 
route is different. As.should be apparent at this point, the decision phase of the 
north(south,east).service event routine must correspond with the service policies and 


priority rules discussed earlier. The program listing is presented in Appendix B. 
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Figure 3. 







Is atleast 
one of jn, js, je 
less than N17 ? 


No 







Schedule 
new.north.unit event routine 
now. 







Reschedule this event routine 
foranext arrival. 





END 


Flow chart for north.unit.generator event routine: Jn 1s the number of 
north.units that have passed the intersection. js 1s for south.units and je 
1s for east.units. Nl is the minimum desired number of units to be sim- 
ulated from each direction and 1s selected by a user. This routine will be 
rescheduled as long as at least one of jn, js or je 1s less than NI and will 
stop being rescheduled if all three values are greater than or equal to 
NI. 


Create 
a north unit now. 








Yes Is any unit served 
or waiting in the 


north queue ? 


No 
Yes No 
Next.north.route = 1. Next.north.route = 2. 


Is there Is there 


blocking to blocking to 
north route! ? north route2 ? 










Yes Yes 





No 

North.run = 1 and ee North.run =2 and 

next.north.route =0. j ; next.north.route = 0. ; 

routine. { 

i 

| 

‘ File this unitin the 
Return i 


Figure 4. Flow chart for new.north.unit event routine: U 1s a Uniform (0.1) ran- 
dom variate and is compared with P, for route selection of the unit where 
P, is the probability of taking route]. An arriving unit which 1s created 
in this routine can be served immediately or filed in the queue depending 


upon queue condition, route selection and the other movements at the 


intersection. 


Remove the unit 
from the system. 






Me Is there any unit 
leftin the 
north queue ? 
Yes 
Yes No 
Make Make | 
next.north.route = 1. next.north.route=2 | 
Yes | No 


Make 
north.run = 0 









Is 
next.east.route = 1 
and is that route 
clear ? 






Is 


next.east.route 
=1? 







No No 


Yes 


Figure 5. Flow chart for north.service event routine: U is a Uniform (0,1) random 


variate and P, is the probability of taking routel. In this routine, prior- 
ities for service after a service completion on north route! or 2 are given 
to appropriate routes according to the service policies discussed in Table 
2 through Table 7 in Chapter II. 


pay 









East.run=1land 
schedule 
east.service even: 
routine. 





East.run =1 and 
schedule 
east.service event 
routine. 





Is 


next.south.route 
=1? 

















No 


Yes 





South.run = 1 and 
schedule 
south.service event 
routine. 












Is 
next.south.route 
=e 


No 







South.run = 2 and 
schedule 
south.service event 
routine. 


Yes 











Is 


next.north.route = 2 
and isthat route 


clear ? 














Is 
next.north.route = 1 
and is that route 
clear ? 





No 


North.run =2 and 
schedule 





North.run = 1 and 
schedule 
north.service event 
routine. 






north service event 
routine. 





Figure 6. Flow chart for north.service event routine (continued): The priorities for 
next service after service completion in the south or east are decided in 
the south(east).service routines. Note that in the program the value of 
next.north(south,east).route — is changed to 0 whenever a 


north(south.east).service event routine is scheduled. 


Begin 


Preamble 


Call read.data routine 


Make a design matrix and 
run a simulation under each condition. 


Call initialize routine 


le 


Start asimulation 
Simultaneously simulation 
clock Starts to run. 


Interactions between routunes. 





} Call single run, sum report. | } Call single run, sum report. | run, sum report. 


Compute and print 
SIMUIATION STATISTICS. 


No 


Reset all counters 
Release all arrays 


Return ) 


Figure 7. A macro flow chart for the model: N is the smallest number of 
north(south,east) units that have passed the intersection and NI (se- 
lected by the user) is the desired number of units to be simulated. R and 
NR indicate the number of replications completed and the desired num- 
ber of replications respectively. C and NC are used to check if all simu- 
lation conditions are examined where C indicates the number of 
conditions that have been examined, and NC is the total number of 


conditions to be examined. See Table § for the simulation conditions. 
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D. RESULTS 
1. Statistical Considerations 

When the combinations of many different input specifications to a simulation 
are to be evaluated, as in this problem, we must first decide upon which combinations 
of factors we wish to spend our limited simulation resources. So, manv preliminary ex- 
periments using various combinations of input parameters were investigated, and those 
conditions under which the system was not stable were used as boundary conditions in 
the experimental design. As a result, a finite number of discrete level experiments were 
conducted as a series of essentially independent tests in looking for trends in each di- 
rection. We used 100 replications within each run to achieve sufficient precision for each 
response estimate. 

In order to reduce the variance among simulation runs, common random num- 
ber sequences were used for all the simulations of different factor combinations, [Ref. 
11: pp.250-252]. Three separate sequences were used for each simulation (interarrival 
time, service time, and route selection). The random number sequences were continued 
to the next replication, so that there was no overlap of the random number sequences 
and we could therefore assume the replications for each set of factors were independent 
of each other. 

2. Simulation Outputs 

Whenever we have one replication of a simulation, as discussed before, a par- 
ticular realization of the output sequences will be generated. Examples of random se- 
quences of interest would be a running average of successive waiting times and a running 
average of queue length experienced by the first y north(south,east) units for r = 
],...,.N1, where NI is the desired number of units to be simulated from each direction. 
NI is a user input and must be large enough to provide stability. 

For this study, the following sequences of interest were collected and computed 


in the single.run.report and sum.report routines. 


W., : Waiting time for the r” departing unit from the /” direction, where / can have 
values north, south or east. 


T,, : Departure time for the r® unit from the 7” direction. 


AW., : The running average of the first yr waiting times obtained by the formula: 


The number of units waiting in the /” direction at the r” unit’s departure. 
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a 


AW,,=— ) W, 


if (2 


p—Woniniesontn.east and r= l,....N1. (3418) 


v4 


AQ,, : The running average of the number of units waiting in the / direction during 


the time 0<7r</T,,. It is a function of time and a random sequence, but it can also be 
viewed as a sample of random processes defined on continuous time or as a function or 
integral of such a process. For example, if we let Q(2) be the number of units waiting in 
the i* direction’s queue at time 7, then {Q(1); > 0} 1s a random process with a contin- 
uous parameter ¢. For :>0, AQ,, can be expressed by the following formula, [Ref: 11: 
Be) 2-273): 


Ti, 
AQ, = |, O(u)du, i=north,south,east andr =1,...,N1. (3.2) 
ir 


which is easily computed by the ACCUMULATE function in Simscript II.5, [Ref. 7: 
Be) 2-273]. 

The sequences of V,,, AW,, ,AQ,, and the process Q(t) for the north direction 
and 100 units of N1 with skewed interarrival and service time distributions, traffic 
intensity=0.5 and P, =90.75 are shown in Figure 8. Note that this set of simulation 
conditions 1s also used for further system output displays. 

3. Characteristics of Output Sequences 

As we discussed above, the outputs of a simulation are sequences of random 
variables. In this study we are interested in certain probability distribution characteristics 
and our objective is to estimate them via this simulation since they are unknown and can 
not be computed with theoretical analysis. 

Assume we have Af independent observations of the sequence {AlV,,} and 
{AQ,,}, /=north,south,east and r= 1,...,N1 through independent replications of a sim- 
ulation. To obtain independent replications of a simulation, we generate independent 
sequences of random variables to drive the simulation. These sequences are recursively 
generated, as discussed previously in this section. 

To estimate the mean values and obtain confidence intervals of the observations 
from sets of sample sequences (such as depicted in Figures 9 and 10) we let u aw and 


oAW,, be the mean and variance sequences of AW,, and H4o,, and of be those of AQ,,, 


0; 
IKe:, 
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Figure 8. Output sequences obtained from a single simulation run 


jie = LA '= north,south,east and r= 1,...,NI. (333) 
Hag, = EVAGw; i= north-south east. and = laa (3.4) 
a", w= EUAM,, — 4, w,dd i= north,south,east and r= 1,...,N]. (3.5) 
2 2 
ao. =E[(AQ,,-—pu,, )], i= north,south,east and r= 1,...,NI1. BES) 
AQ, : AO e 
A — l 


A Wie = 


M 
a » A Winm,  ¢=north,south,east and r= 1,...,NI. (3.7) 
nil 


PAW, 


Thus, the sequence of standard point estimators of w,,, is the sequence of 


sample means 1, ryt (3.7) where M is the desired number of replications of the simu- 


oe 


lations and AWV,,,, is the value of the m” replication of AW. A confidence interval for 


the sequence y,,, is provided with the t-statistic as follows, [Ref. Ll: pp.282-283]. 


A 


AM 
: : | Te 
Sim," Sau, D4 Wim —- AW, J, (3.8) 





i= north,south,east and r = 1,...Nl and M = number of replications 





Sy v= c, Ww. i= north,south,east and r = I,...,NI1. (329) 
A S. Wi, : 
Vary , y= va 2 = Wer, southeast and r = l....,N1. (3.10) 


Given the above definitions a random variable, 


A 
BA Wis 7 P 4 Wie 
Sy W’. 


¢,? 


Mir 





has approximately a t-distribution with ‘’f— 1 degrees of freedom and we have 


A 


Law —H 
Prob\ty,_(a/2) < ar: Ain Lago be OG eae — OL, (S221) 
AW, 


iF 


Wee 





where /=north,south,east, r=1,...N1, M= number of replications, and t,_,(a/2) is the 
o/2 quantile of the t-distribution with 4 — | degrees of freedom. Since the t-distribution 


is symmetric about zero, we obtain, 








A Si Mie A = ue ft a 
a — nN Oldhesouth ecastuandure—lewsN 1. 
A S, Wie A AW,, 
Therefore the two values, y, Ww. ty_\(1 — a/2) Vir and yu, w,, + ty—\(1 — o/2) yi are 


100(1 - «)% confidence intervals for u,,,. Given the previously discussed paraineter 
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conditions, the realizations of 4,,, and the confidence intervals of u,,, obtained from 
the 6 sequences in Figure 9 and Figure 10 are depicted in Figure 11 and 12 for the north 
and south directions respectively. 
A hy . A : e 

Let Lyy,.., b€ the 500” value of the sequence yu aw,,) the mean estimator of the 
running average of successive waiting times experienced by the first 500 units and let 
FAW, soo be the standard deviation estimator of u AW, son” The values of p AW, con’ PAW, soe ob- 
tained from various simulation conditions are shown in Tables 16 through 18 of Ap- 
pendix A. The values of p AW, eq) 108 different parameter sets and different directions will 


be analyzed in detail in Chapter V. 
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300 400 500 





Figure 9. Six realizations of an output sequence (north queue): This figure shows 
six realizations of the sequence {AW,,} for value 1 corresponding to the 
north direction. We can see the random appearance and the variation 
from sequence to sequence. The distribution characteristics of the ran- 
dom sequences such as these are what we are interested in estimating. 
These characteristics are estimated from the sequences in the figure and 


are Shown 1n Figure 11. 





Figure 10. Six realizations of an output sequence (south queue): This figure shows 


six realizations of the sequence {AW ,,} for value 1 corresponding to the 
South direction. The sequences in this figure also show random ap- 
pearance and variation from sequence to sequence. The distribution 
characteristics of the random sequences are estimated and shown in 


reure ))2: 


4. Stability and Convergence Considerations. 

Note that there are three separate queues in the system. As we see in Figure 
11 and Figure 12, the sequence (Haw, seems to CONVETgE tO Hay, in the north and south 
queues under the given traffic conditions. We found out that the east queue is not stable 
under the same conditions. We Know that it makes a difference in simulation difficulty 
if the sequence of interest settles down and converges to a stationary distribution after 
an initial transient phase. If the sequence settles down, then we may only need to con- 
sider the distribution that exists after the transient phase and not a separate distribution 
for each value in the sequence, [Ref. 6: p.II].22]. 

To test convergence to the stationary distribution, we made preliminary, inde- 


pendent replications of the simulation over a range of r (number of units simulated) 
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Figure tl. Mean and C.I sequences of waiting time (north queue): This figure 
shows a realization of {nu Ai) and confidence intervals of {nu Ae) for 
value 1 corresponding to the North direction. The sequence (Figy ) 1S 
a random sequence obtained from 6 sequences of AW’,, in Figure oF 
Since EAE op, =H,y, » it fluctuates about the converging sequence of 
means. [he larger the value of M the narrower the fluctuations and the 
closer the approximation to the sequence eo Je Hence, bv taking a 
large value M (100 in this thesis) a reasonable estimate can be made 
on a convergence point for the sequence erg 3. Notice that from this 
figure a reasonable judgement can be made as to the extent of the 
transient phase, whereas this would be impossible from any individual 
output sequence. In this figure a reasonable estimate for the end of the 


transient phase appears to be about the first 100 waiting times. 


sufficient to observe convergence. The sequence {u,,,, } was estimated by generating a 
sequence of sample means as discussed before. The most straightforward convergence 
test procedure would be to collect and plot the whole sequence and check the conver- 


gence by examining the trends. In practice, however, such an investigation requires more 
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Figure 12. Mean and C.I sequences of waiting time (south queue): This figure 
shows a realization of ae and confidence intervals of ce for 
value 1 corresponding to the south direction. The sequence {ligy } is a 
random sequence obtained from 6 sequences of AlV,, in Figure 10. In 
this figure we can see that the south queue is stable and a reasonable 
estimate for the end of the transient phase appears to be about the first 


300 waiting times. 


. ; hi 
computing time and effort than would be reasonable. Therefore we collected every 10 


data point after the first 300 points of a sample mean sequence, such as the one seen in 
Figure 11, to save our limited simulation resources. To examine convergence of the re- 
sulting sequence we manually compared the values over the sequence range. 

The convergence tests revealed the system is partially stable in the north and 
south queues for Model | under the condition that traffic intensity=0.5, P, = 0.75 and 
the ratio of mean interarrival time for east over mean interarrival time for north or south 
is 1. To illustrate the east direction’s unstable queue, we collected 6 complete sequences 
of AV,,,,, and plotted them in Figure 13. Note the upward trend in the sequences. As 


we increase the ratio; (mean interarrival time for east) / (mean interarrival time for 
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north), from | to 1.3 by increments of 0.1, the east queue becomes more stable (Figure 
14). 





Figure 13. Six realizations of an output sequence (east quee); This figure shows 
six realizations of the sequence {AW’,,} for value 1 corresponding to the 
east direction. The sequences show an upward trend and indicate there 
is no steady state distribution for the mean of the sequences. The 
characteristics of each data point in the mean sequence can be esti- 


mated from the corresponding data points in the six AW,, sequences. 


5. Data Collection and Analysis Tools 

We have powerful tools for data collection and analysis such as “print” and 
“accumulate/tally” statements in PC-simscript. These statements can be used for data 
collection and some statistical calculations. Other software packages such as 
STATGRAF in PC can be used for more sophisticated data analysis. In addition to 
these tools for use on the personal computer, we used IBM 3033AP mainframe support 
to efficiently manipulate much of the data. Data was collected on floppy diskettes and 
transferred to the IBM mainframe with the 3270 PC’s Host-Supported File Transfer 


program. Once transferred to the mainframe, we edited the data file prior to analyzing 
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the data with GRAFSTAT. GRAFSTAT graphical displays were extremely useful for 
checking the characteristics of large data sets such as required for Figure 9 and Figures 
16 through 22 in Chapter V. 
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Figure 14. Stability test on the east queue: We tested the stability of the east 
queue by increasing the ratio of mean interarrival time for the east over 
the mean interarrival time for the north from 1.0 to 1.3. The estimators 
of ELAIV,,] = (Hy) is 2 the sequence of the sample means of 
{Al}V,,} for value 1 corresponding to the east direction. The east queue 


becomes more stable as the “ratio” increases. 
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IV. SIMULATION MODEL 2 


A. ASSUMPTIONS AND CHANGES TO THE MODEL 

After developing a model of the system using the described assumptions, we propose 
and simulate a road structure to reduce the restrictions of the current road system. 
Figure 15 depicts the proposed system which will be addressed as Model 2. It 1s obvious 
that Model 2 allows more combinations of routes, so that traffic performance at the 
intersection should be increased. Actually 4 combinations of routes, (0,1,1), (2,0,1), 
(2,1,0), and (2,1,1), are added to the 14 cases of Model 1. In case of blockages, the pni- 
orities for next service are shown in Table 9. Note that in Model 2, each route only 
blocks the route that merges with it in a common direction. All other assumptions of 


Model 2 remain the same as in Model 1. 


Table 9. PRIORITIES AT THE END OF SERVICE OF A BLOCKING ROUTE 
(MODEL 2): See Table 1 and explanation of Table arrangement pro- 


vided on page 7 for comparison with Model 1. 


Route which is about ae 
. Next Priority, in order, at the end of blocking 
to finish service 


north routel east routel, (north routel, north route2) 
north route2 south route2, (north routel, north route2) 


south routel east route2, (south routel, south route2) 


south route2 north route2, (south routel, south route2) 


east routel north routel, (east routel, east route2) 


east route2 south routel, (east routel, east route2) 
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Figure 15. 


East route2 East route! 


North 
route2 


North route1 


South route2 


- South routel 


-< 


A proposed modification of the current system (Model 2): Here several 
blockages have been removed. In particular traffic coming from east 
routel is never blocked by traffic coming from south routel though it 
sull has to merge with traffic from north routel. Similarly, traffic on 
north route2 is no longer blocked bv either traffic on east routel or 
south routel, though it still can be blocked by traffic from South 
route2. Finally, traffic on south routel is no longer blocked by traffic 
on north route2 or east routel but can still be blocked by traffic from 


east route., 


B. INPUTS TO THE MODEL 


Inputs are the same as in Model 1. 


C. MODEL IMPLEMENTATION 

Since the only difference between Model 1 and Model 2 are the additional route 
combinations, all other events and routines remain the same as the first model except the 
events related to evoking service. These routines are new.north(south,east).unit and 
north(south,east).service event routines. As an example, consider the route combination 
(2,1,0); it is not allowed in Model 1 because of the movements of north route2 and south 
route | are in conflict, however, the two movements are allowed in Model 2. Thus the 
decision phase of the event routines are revised for the new model as required to allow 


the additional route combinations. The program listing is shown in Appendix C. 


D. MODEL VERIFICATION AND VALIDATION 

Having set up the service policies, system assumptions and developed the simulation 
model, the next step is to verify the model; that is, to detect and correct discrepancies 
between the intended and actual execution processing performed by the model. 

A simulation model of a complicated process such as this model can not be de- 
bugged or verified unless comprehensive diagnostic outputs, which completely describe 
the operation of the simulated svstem, are provided. Diagnostics are needed only for 
debugging, being much too cumbersome and costly to use in full-scale runs. The most 
powerful and economic diagnostic tool used in this model is output which describes each 
simulated event on a separate line with almost all relevant information about each unit's 
status. The data collected in the single.run.report routine is also used for verification 
of the model. The sample sections of this diagnostic and the data collected in the 
single.run.report routine are shown in Tables 13 through 15 of Appendix A. Another 
verification tool used was to examine the logic of the flow charts discussed earlier. 

To validate the model as a sufficiently accurate representation of the intersection 
operations for the purpose of the study, that is, to answer the question of whether there 
is a correspondence between the model and the physical intersection system, would re- 
quire the models results be compared with that of the real system in a given situation. 
It would be impossible to create the physical system conditions simulated by the model 
because of the many restricted military situations in R.O.K. We can not exercise the 
movements of all forces to simulate wartime traffic given the current situation. In spite 
of the restrictions, manual checks of the model flow charts and careful examination of 


the assumptions and service policies may be conducted to confirm that the model cor- 
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responds to the actual system. If we assume that military police control the traffic at the 
three way intersection, the general assumptions and service policies in the models do not 
correspond to the real world, however, that 1s not a problem since we are interested in 
overall traffic performance at the intersection and service priority issues do not affect 
our study objective. 
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V. ANALYSIS 


A. SIMULATION RESULTS 
1. Current Road System (Model 1) 

We present a graphical approach in organizing and presenting the simulation 
results, [Ref. 6: p. VIII.3-23]. The interest 1s showing, in a simple graphical form, how 
the chosen quantification, 1.e. mean waiting times at the intersection, vary with factors 
such as traffic intensities, service(unit) length and interarrival time distributions. Figures 
16 through 18 show the results of the first model with factor P, = 0.75, 0.50 and 0.25 
respectively. 

In Figure 16, we have six sub-figures; three in the upper half with mean service 
times of 2.5 for all directions and three in the lower half with mean service times of 5.0 
for all directions. Each cell in the sub-figure is indexed by the traffic intensity and the 
combinations of interarrival and service time distributions defined here to be SS, SR, 
RS, and RR, (see Figure 16 for an explanation of these terms) and contains the esti- 
mated value of yu, W, e00 (see Equations 3.7 and 3.9). In Figures 17 and 18, we can observe 
mean waiting times of the three directions for all combinations of interarrival and service 
time distributions, given traffic intensitv=0.5 and mean service time= 2.5. In addition, 
Figure 22 shows the changes to the mean waiting times of the east queue together with 
changes to the mean interarrival time ratios for the three directions, given Model |] with 
P, = 0.75, traffic intensity = 0.5, and mean service time= 2.5. Note that the waiting times 
in Figures 17 and 18 should be compared with the waiting times on the graphs 1n Figure 
16 with mean service time 2.5 and traffic intensity of 0.5. Note too that the vertical 
scales in Figures 16 through 18 are different. It would be preferable to use the same scale 
to give a depiction of the relative difference of the mean waiting time as simulation pa- 
rameters change, however, this would require either a very large display or a transfor- 
mation of the data (such as a log scaling). Neither of these were desirable, thus we need 
to look at the vertical scales when comparing the estimators between sub-figures. 

Hence, it 1s relatively easy to visually inspect Figure 16 to determine the effect 
of the mean service time (by looking at the figure’s columns), the effects of traffic in- 
tensity (bv comparing the rows in a sub-figure), and the effect of the combinations of 


interarrival and service time distributions (by looking at the columns in a sub-figure). 
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However, in order to observe the effects of changes to the probability that the first unit 
in the queue selects routel, P,, we need Figures 16 through 18 side-by-side. 

It is immediately clear in Figures 16 through 18 that traffic intensity p has a 
major effect on the mean waiting time and, in fact, it is known from theory that the 
mean waiting time in many situations 1s proportional to the reciprocal of (1-p), [Ref. 6: 
p.VIII.8]. Beyond this, it appears to be generally true that the regular-regular combi- 
nation of interarrival and service times at the right most position 1n the six skyscraper 
plots of Figure 16 is best for traffic performance; on the other hand, the skewed-skewed 
combination at the left most position in the six plots is the worst case for all values of 
traffic intensity. In addition, as the mean service time doubles the mean waiting time 
also doubles. Finally, the effects of increasing mean service time from 2.5 to 5.0 together 
with other factor combinations turns out to be consistent with different values of P,. 

To analyze the effects of the factor P, requires comparing Figure 16 (P,=90.75) 
through Figure 18 (P,=9.25). For instance, the mean waiting time of the east direction 
is almost five times that of the other directions in Figure 16 with P,=90.75. On the other 
hand, Figure 17 with P, = 0.5 shows that the mean waiting times for the south direction 
are relatively higher than the north and east directions. This is because the frequency 
of movements on each route varies as the value of P, changes in conjunction with dif- 
ferent route movement conditions. 

Using Figure 22 for P, = 1 and traffic intensity = 0.5, we found that the system 


is stable in the north and south directions but can be unstable in the east direction if 


mean interarrival time of east units 
mean interarrival time of north or south units ~ 


— 
2 


A more sophisticated analysis of the results, for instance the test to see if there 
still remains a difference in mean waiting times of three directions due to P, after ad- 
justing other factors, 1s not done in this paper, [Ref. 6: p. VIII.9]. 

2. Proposed Road System (Model 2) 

The results of the simulation for Model 2 are shown in Figures 19 through 21 
and the waiting times in these figures are compared with the corresponding values of 
Figures 16 through 18. The effects of the factors found in the first model are consistent 
in this model. We can expect that the overall traffic performance in Model 2 would in- 
crease because the new system allows more combinations of movements, that 1s, the re- 
Strictions to movements are reduced. The benefits of the new svstem for the north 


direction removes two possible blocking routes, east route! and south routel, which are 
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blocking units in the first model for north route2. As for the south direction, east route] 
and north route2 no longer block south routel. Finally in the case of the east direction, 
east route! is no longer blocked by south routel or north route2. Thus, overall move- 
ment conditions are improved and waiting times, as would be expected, are decreased for 


the three directions in this model. 


B. COMPARISON OF THE TWO MODELS 

In order to compare the results of the two simulations, the sample means and 
standard deviations of the first 500 average (running) waiting times for each direction 
(see equations 3.7 and 3.9) are tabulated in Tables 19 through 21 of Appendix A, so that 
a direct comparison of the waiting times for the two models is possible on a pair-wise 
basis. The fact that the new model 1s more efficient than the current model in traffic 
performance is evident. However, the criteria to prefer the new system to the old could 
be more than just the relative differences in time. Numbers of additional units deploved 
per unit of time, the decrease in forces susceptible to enemy air and ground ambushes, 


are a few of the criteria that could be used to compare the two systems. 
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Model 1 and P, = 0.75 for all directions 
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Figure 16. Effects of major factors for Model I: Each bar in the cell shows the 
mean estimate of the first 500 (N1=500) waiting times for each di- 
rections queue. The notation SS, SR, RS, and RR on the X axis rep- 
resent combinations of interarrival and service time distributions where 
S and R means ‘skewed’ and ‘regular’ interarrival or service times re- 
spectively. For instance, the notation SR indicates ‘skewed’ interarrival 
and ‘regular’ service time distributions were used for the simulation. 
These are achieved by using different Gamma shape parameters 4,, k, 
for the interarrival and service times. The effects of the traffic intensity 
and the combination of interarrival and service time distributions are 
obvious and consistent in the six skyscraper plots. The mean waiting 
times in the east direction are very large and unstable while the other 
two directions are stable. This is discussed further in Figure 22. Note 
that the scales are different between the six sub-figures. 


47 


Mode! 1 and P, = 0.50 for all directions 
Traffic Intensity = 0.5 
Meon Service Time = 2.5 
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Figure 17. Effects of probability of taking route1(0.50) for Model 1: See Figure 
16 for the explanation of SS, SR, RS, and RR in the plot. The mean 
Waiting times in this plot correspond to a mean Service time 2.5 and 
traffic intensitv of 0.5 for each direction. The effects of decreasing P, 
from 0.75 in Figure 16 to 0.50 in this figure 1s to decrease the east’s 
Waiting time and to increase the south’s waiting time. The decrease in 
the east’s waiting time is because the probability that east units take 
east route2 increases to 0.50 from 0.25 in Figure 16 and it 1s easier for 
east units to move to east route2 than east routel since east routel has 
three blocking routes. On the other hand, there 1s only one blocking 
route for east route2. The reason for the increase in south’s waiting 
time 1s because more north units select north route2 and more east 
units select east route2 which combine to increase the difficulty for 
south units to move under P,=0.50 vs P, =0.75. 
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Model 1 and P, = 0.25 for all directions 
Traffic Intensity = 0.5 
Mean Service Time = 2.5 
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Figure 18. Effects of probability of taking route1(0.25) for Model I: See Figure 
16 for the explanation of SS, SR, RS, and RR in the plot. The mean 
waiting times in this figure are compared with the corresponding values 
in the third row of the upper three plots of Figure 16 and against Figure 
17. The major effects of decreasing the value of P, from 0.50 in Figure 
17 to 0.25 in this figure is to decrease east’s and south’s waiting times. 
The decrease in east’s mean waiting time 1s because now east units se- 
lect east route 2 with probability 0.75 and east route2 1s only blocked 
by south routel, which 1s selected by south units with probability 0.25. 
That is, the probability that an arriving east unit is blocked becomes 
lower than in Figures 16 or 17. The decrease in the south waiting times 
can be contributed to the difference in the probability of taking each 
route and the interactions of movements on the routes. 
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Medel 2 and P, = Q.75 for all directions 


Mean Service Time = 2.5 


*  Woiting Time 
Os : 12 3 4 
“) 

Woiting Time 


Oe) 


<r, 
rival Service Time 
pistribution 


—— RR ER 

sR OS. 6 se OS 
enval,seric® Tim ’ crarrival,Service TES 
" Gistribution 


5S 
intete Distribution 


x 
SS 


ss 
interor 


Mean Service Time = 5.0 
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Figure 19. Effects of major factors for Model 2: See Figure 16 for the explanation 
of SS, SR, RS, and RR in the plots. The waiting times on the north, 
south and east directions decreased almost 2, 6 and more than 10 times 
respectively from the corresponding waiting times of Figure 16 in the 
first model. The differences in the rates are because the performance 
benefits of the new svstem are different for each direction. For example 
with P, = 0.75, the benefits of new east routel and north route2, are 
different since approximately 75% of the east units use east routel and 
only about 25% of the north units use north route2. However, the 
overall effects of traffic intensity and mean service time remain the 


same as 1n Model 1. 
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Figure 20. 


Model 2 and P, = 0.50 for all directions 
Traffic Intensity = 0.5 
Mean Service Time = 2.5 
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Effects of probability of taking route1(0.50) for Model 2: See Figure 
16 for the explanation of SS, SR, RS, and RR in the plot. The mean 
Waiting times in this figure can be compared with the corresponding 
values in Figure 17 to see Model 2’s effect. The corresponding values 
in the third row of each direction for the upper three plots in Figure 19 
depict the effects of changes to P, in Model 2. Since each direction has 
only one blocking route and P, = 0.5, it is reasonable to expect that the 
mean waiting times of the three directions are approximatelv equal. 
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Figure 21. 


Model 2 and P, = 0.25 for all directions 
Traffic Intensity = 0.5 
Mean Service Time = 2.5 
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Effects of probability of taking route1(0.25) for Model 2: See Figure 
16 for the explanation of SS, SR, RS, and RR in the plot. It 1s again 
obvious that the mean waiting times decrease in comparison with the 
corresponding case of the first model (Figure 18). The south and north 
have longer waiting times than the east because north and south units 
select north route2 and south route2 with probability 0.75. On the 
other hand, east units interact only with the north and south units 
which take routes blocking east units with probability 0.25. 
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Model 1 and P, = 0.75 for all directions 


Traffic Intensity = 0.5 
Mean Service Time = 2.5 
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Figure 22. Stability changes on east queue (Model 1): This figure shows the 
specific characteristics of the intersection given Model 1 with P,=0.75 
and traffic intensity = 0.5. For these given input parameters, this figure 
shows changes of stability on the east queue given changes to the mean 
interarrival time ratios for the three directions. The changes are con- 
sistent with any combination of interarrival and service time distrib- 
utions, 1.e. SS, SR, RS, and RR (see Figure 16 for an explanation of 
So oes >and Rie The results in this figure are from skewed 
interarrival and service time distributions (see page 21). The numbers 
on the X axis; 1.0, 1.1... represent the ratio of (mean interarrival time 
of east units / mean interarrival time of north or south units). As seen 
in Figure 16, the east queue is unstable if the arrival rate from the east 
is the same as that of north and south, and becomes stable as the ar- 
rival rate from the east decreases with respect to the north or south 
arrival rate. North’s and south’s stability was not affected by the ratio 
changes. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


This research has demonstrated the effects of traffic conditions on traffic flow at two 
kinds of three way intersections. We have seen how traffic along the roads of interest 
varies as the conditions change. The question of how we interpret mathematical terms 
such as Gamma distributions and other random variables in the real world still remains 
uncertain. That is, how do we know actual interarrival and service times would be dis- 
tributed with certain characteristics? For example, we have assumed a Gamma distrib- 
ution for interarrival and service times in this model, however this assumption may be 
completely different from the real world situation. 

As for service time, it seems as if many factors affect service in wartime in addition 
to the physical restrictions of a road system. As mentioned before, major problems such 
as refugee flow along the roads should be studied separately. But we can assume a cer- 
tain distribution for service time for the purpose of this study. 

For interarrival times, however, we can collect the appropriate data by conducting 
realistic war game simulation. As the war game situation changes, the target level forces 
are assumed to complete their movement according to the combat model’s movement 
schedule, (Ref. 12: p. 3.15-23]. The movement speed is assumed to vary according to the 
characteristics of the theater. However, the current models do not consider the road 
movement interaction of forces which may create critical delavs in combat force oper- 
ations. Therefore we must consider the road movement interaction of combat forces 
given Various Wartime conditions and find out whether the movements of the target level 
forces in the wargame simulation can be completed in a timely manner in the real world. 
This study enables us to analyze the delavs in combat force operations by considering 
the interaction of convoys which have to traverse an intersection of interest. If we create 
a road march plan for the forces according to a war game simulation situation then we 
can collect the expected arrival times of each unit at the intersection of interest. The 
expected arrival times at an intersection of interest obtained from this procedure could 
be developed into a distribution for interarrival times. 

Thus we can apply this model to the real world with the given assumptions and can 
determine expected traffic flows in the area of interest. In this thesis, we study only one 
three way intersection in an attempt to describe its overall traffic performance. It would 
be especially useful to apply this model to a three way intersection which is believed to 


be a traffic choke point in a specific area of interest. Further studies about more ag- 
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gregated road systems along with war game simulation are needed to further understand 
the potential for wartime traffic problems in the Korean theater. 
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APPENDIX A. SAMPLE OUTPUT OF THE SIMULATION MODELS. 


Table 10. SAMPLE SECTION OF DIAGNOSTIC PRINTOUT OF MODEL 
1: Simulation conditions; P, =0.75, mean service time = 2.5, Tl(traffic 
intensity) =0.5, interarrival and service time distribution=SS (see Table 
8 and Figure 16). 


NNU: new. north. unit, NSU: new. south. unit, NEU: new. east. unit 
NS: north. service, SS: south. service, ES: east. service 
U: a uniform random variable used for the next route selection. 


Movements: movements at the intersection after the left event. 
N,S,E: north, south, and east direction. 

Queue: number of units in each queue after the left event. 
Next route: routes selected for the first unit in each queue. 
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Table 11. SAMPLE SECTION OF SINGLE.RUN.REPORT OF MODEL 1. 


arrival time of a unit. 

queue length after an arrival. 

queue length before departure of a unit. 
waiting time of each unit. 

departure time of a unit. 


time of service completion. 
route number used by a unit. 
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Table 12. SAMPLE SECTION OF SUM.REPORT OF MODEL 1. 






AQn: average queue length(running) at each north unit's departure. 
AWn: average waiting time(running) at each north unit's departure. 


AQs and AWs are for south units, AQe and AWe are for east units. 


No AQn AWn AQs AWs AQe AWe 
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4 . 89 2.79 .44 Salis Ya) ih 9.26 
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Table 13. SAMPLE SECTION OF DIAGNOSTIC PRINTOUT OF MODEL 
2.: Simulation conditions; P, = 0.75, mean service time = 2.5, Tl(traffic 
intensity) = 0.5, interarrival and service time distribution=SS (see Table 
8 and Figure 16). 


NNU: new. north. unit, NSU: new. south. unit, NEU: new. east. unit 
NS: north. service, SS: south. service, ES: east. service 
U: a uniform random variable used for the next route selection. 


Movements: movements at the intersection after the left event. 
N,S,E: north, south, and east direction. 

Queue: number of units in each queue after the left event. 
Next route: routes selected for the first unit in each queue. 





Movements Queue Next route 
Event Time U (Nemo £) CN Sa E) CNS) 
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Table 14. SAMPLE SECTION OF SINGLE.RUN.REPORT OF MODEL 2. 


arrival time of a unit. 
queue length after an arrival. 
queue length before departure of a unit. 


waiting time of each unit. 
departure time of a unit. 
time of service completion. 
route number used by a unit. 
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Table 15. SAMPLE SECTION OF SUM.REPORT OF MODEL 2. 


AQn: average queue length(running) at each north unit's departure. 


AWn: average waiting time(running) at each north unit's departure. 
AQs and AWs are for south units, AQe and AWe are for east units. 
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Table 16. SIMULATION OUTPUTS FOR VARIOUS CONDITIONS OF 
MODEL |! 










probability taking routel. 
MS: mean service time. 

TI: traffic intensity(see Table 8). 

A.S: interarrival and service time distribution (see Figure 16). 
WT: mean waiting time (see equation (3.7) in chapter III. ) 
std.dev. of WI(see equation (3.8) in chapter III. ) 





2.39 Wee 
1.82 W20. 
10. 9.0 265)  2Greae 3.21 159.92 Tomes 
-96 320: 









62 


Table 17. CONTINUED FROM THE PREVIOUS TABLE. 
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Table 18. CONTINUED FROM THE PREVIOUS TABLE. 
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Table 19. SIMULATION OUTPUTS FOR VARIOUS CONDITIONS OF 
MODEL 2 


probability taking routel. 

mean service time. 

traffic intensity(see Table 8). 

interarrival and service time distribution (see Figure 16). 
mean waiting time (see equation (3.7) in chapter III.) 
std.dev. of WTI(see equation (3.8) in chapter III.) 
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Table 20. CONTINUED FROM THE PREVIOUS TABLE. 
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Table 21. CONTINUED FROM THE PREVIOUS TABLE. 
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APPENDIX B. A PROGRAM LISTING OF MODEL 1. 


PREAMBLE ''an intersection simulation 

normally, mode is undefined 

event notices include north. unit. generator, 
south. unit. generator, 
east. unit. generator, 


new. north. unit, 
new. south. unit, 
new. east. unit, 


north. service, 
south. service, 
east. service 


every north.service has an north. pointer 
every south. service has an south. pointer 
every east.service has an east. pointer 


define north. pointer, 
south. pointer, | 
east. pointer as pointer variables 


temporary entities 

every unit may belong to the north. queue 
may belong to the south. queue 
may belong to the east. queue 


the system owns the north. queue , 
the south. queue, 
the east. queue 


define north. queue, 
south. queue 
and east. queue as fifo sets 


define north. run, 
south. run, 
east. ein, 
next. north. route, 
next. south. route, 
Next. Cast. Toure. 
Nl, 
replication,NR, 
in, jn,kn,is,js,ks,ie,je,ke, 
saveseedl ,saveseed2,saveseed3 as integer variables 


define AQ.north.nl, AQ.south.nl, AQ. east.nl, 
AW.north.nl, AW. south.nl, AW. east.nl, 
alphal, betal, alpha2, beta2, 
Pl, current. time, 
Sl,s2,s3,rl,r2,r3 as real variables 
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define 


define 


tally 
tally 


tally 
tally 


tally 
tally 
accumu 
accumu 


accumu 


priori 


END =p 


MAIN '' 
define 


NOrei di reseuen. dir,ecast. dir 
as l-dimensional integer array 


arrival. time. of. north. unit, 
departure. time. of. north. unit, 
end. time. of. north. service, 
Waoleine, bime, Of north, unit, 
north. queue. lengthl, 

north. queue. length2, 


ameivei tame, of. south. unit, 
departure. time. of. south. unit, 
end. time. of. south. service, 
waiting. time. of. south. unit, 
south. queue. lengthl, 

south. queue. length2, 


arrival. time. of. east. unit, 
departure. time. of. east. unit, 
end. time. of. east. service, 
waiting. time. of. east. unit, 
east. queue. lengthl, 

and east. queue. length2, 

AQ. of. north. unit, AQ.of. south. unit, AQ. of. east. unit, 
AW. of.north. unit, AW.of.south. unit, AW. of. east. unit 
as l1-dimensional real array 


ave. que.north.nl as the mean and 
sd. ave. que.north.nl as the std.dev of AQ. north.nl 
ave.wt.north.nl as the mean and 
sd. ave.wt.north.nl as the std.dev of AW. north. nl 


ave. que.south.nl as the mean and 
sd. ave. que. south.nl as the std.dev of AQ. south. nl 
ave.wt.south.nl as the mean and 
sd. ave.wt.south.nl as the std.dev of AW. south.nl 


ave. que.east.nl as the mean and 
sd. ave. que. east.nl as the std.dev of AQ. east.nl 
ave.wt.east.nl as the mean and 
sd. ave.wt.east.nl as the std.dev of AW. east.nl 


late AQ.north. unit as the mean of n.north. queue 
late AQ. south.unit as the mean of n. south. queue 
late AQ. east.unit as the mean of n. east. queue 
ty order is north. unit. generator, 
south. unit. generator, 
east. unit. generator 
reamble 


program starts execution here. 
S. que. north, 
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7 Ws HOt 

. que. south, 

Wr sOutur 

. que. east, 

.Wt. east, 

inverse,traffic.intensity as real variables 


nNnnNnNN NW 


call read. data 


saveseedl = seed. v(s1) 
saveseed2 = seed. v(s2) 
saveseed3 = seed. v(s3) 
for Pl = 0.25 to 0. 75) byaueze 
do 
also for alpha2 = 2.5 to 15,0372 > 
do 
also for inverse = 2 to 3 by .5 
do 


also fo 


also fo 


for rep 
do 
le 
le 
le 


le 
le 
le 


traffic. intensity = 1/inverse 
alphal = alpha2*“inverse 


rl = 

r2 = 1 

r3 = l 

r betal = 1 to 4 by 3 
do 

r beta2 = 1 to 4 by 3 
do 


print 2 line thus 
Pal al bl a2 2 eee 
print 2 line with Pl,alphal,betal,alpha2, 
beta2,traffic. intensity thus 


S| esc Teves | Seve Tew | ves eae vesie Weve | vJeeie ve. seve 








lication = 1 to NR 














t seed. v(sl) = saveseedl 
t seed. v(s2) = saveseed2 
t seed.v(s3) = saveseed3 


call initialize 

start simulation 

call single. run. report 
call sum. report 


.north.nl = AQ. of. north. unit(N1) 
.south.nl = AQ. of. south. unit(N1) 
.east.nl = AQ.of. east. unit(N1) 

.north.nl = AW. of. north. unit(N1) 
.south.nl = AW. of. south. unit(N1) 
.east.nl = AW.of. east. unit(N1) 

t time.v = 0 


t saveseedl = seed. v(sl) 
t saveseed2 = seed. v(s2) 
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— etl 


let 


let 


saveseed3 = seed. v(s3) 


Bel 
is 
ie 


| 
ooo 


jn 
js 
je 


neu a 
oo © 


kn 
ks 


0 
0 
ke 0 


if replication = NR and NR > 1 


see 





weve aedede 
Tere deaevke 
always 


mmm me a re em ee me eee ee ee ee 


S.que.north = sd. ave. que. north. nl/sqrt.f(NR - 1) 
S. que. south = sd. ave. que. south. nl/sqrt. f(NR - 1) 
S.que.east = sd. ave. que. east. nl/sqrt. f(NR - 1) 


sd. ave.wt.north.nl/sqrt.f(NR - 1) 
sd. ave.wt. south. nl/sqrt. f(NR - 1) 
sd. ave. wt. east.nl/sqrt.f(NR - 1) 


s.wt. north 
.wt. south 
S.wt. east 


“” 


print estimators of mean queue length (running) 
at the Nlth unit's departure and mean waiting time 
(running) of the Nith unit and std.dev of each. 


print 5 lines with ave. que. north.nl, 
S. que. north, 
ave. que. south.nl, 
S. que. south, 
ave. que. east.nl, 
S. que. east, 
aVveawCewMonrcil. a1,s. Wwe. north, 
ave.wt. south.nl,s.wt. south, 
ave.wt. east.nl,s.wt.east thus 


SSSSS=] | 20 =—=——=—_= == _ == 


main routine for information of data. 














Fede | deveve deve | eveve vex’ dedede weve desde tere dedeos 
Wie dedede dene devede dea aedeve ede device dene device 


release arrival. time. of. north. unit, 


departure. time. of. north. unit, 
end. time. of. north. service, 
Wouudie. CiMe. Of momen. Unit, 
north. queue. length2, 
arrival. time. of. south. unit, 
departure. time. of. south. unit, 
end. time. of. south. service, 
waiting. time. of. south. unit, 
south. queue. length2, 
arrival. time. of. east. unit, 
departure. time. of. east. unit, 
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end. time. of. east. service, 
waiting. time. of. east. unit, 
east. queue. length2, 

AOVOE nocen. unt. 

BOVOE souchs Unmet. 

AQ. of. east. unit, 

AW. of. north. unit, 

AW. of. south. unit, 

AW. of. east. unit, 

north. queue. lengthl, 
south. queue. lengthl, 

east. queue. lengthl, 

north. dir,south. dir,east. dir 


loop 


reset the totals of AQ. north.N1,AQ. south. N1,AQ. east.Nl, 
AW. north. N1,AW. south. N1,AW. east. N1 
loop 
END ~ 'MAIN 


ROUTINE to READ. DATA 
print 4 lines thus 


com mm cm ms we ws ee St ee ee ee ee ee 





see MAIN routine to see other factors 
that are not asked here ! 
print 2 lines thus 
enter the seed stream numberl,2,3 
sl, s2 ,s3 (1 - 10) 
read sl,s2,s3 





print 3 lines thus 
enter the desired number of units 
and replications simulated ! 
N1, NR 
read N1,NR 
END '' READ. DATA 


ROUTINE to INITIALIZE 


let seed. v(sl) = saveseedl 
let seed. v(s2) = saveseed2 
let seed. v(s3) = saveseed3 


schedule a north. unit. generator in 

gamma. f(rl*alphal,betal,sl) minutes 
schedule a south. unit. generator in 

gamma. f(r2*alphal,betal,s1) minutes 
schedule a east. unit. generator in 

gamma. f(r3*alphal,betal,sl) minutes 
print 4 lines thus 

Movements Queue Next route 











2. 


Event Time U (Neo 5) (NYS ,E) Nis Bp) 


reserve arrival.time. of. north. unit, 
departure. time. of. north. unit, 
end. time. of. north. service, 
waiting. time. of. north. unit 
and north. queue. length2 as 2000 











reserve north. queue. lengthl, 
south. queue. lengthl, 
east. queue. lengthl as 2000 


reserve arrival.time. of. south. unit, 
departure. time. of. south. unit, 
end. time. of. south. service, 
waiting. time. of. south. unit 
and south. queue. length2 as 2000 


reserve arrival. time. of. east. unit, 
departure. time. of. east. unit, 
end. time. of. east. service, 
waiting. time. of. east. unit 
and east. queue. length2 as 2000 


reserve north. dir,south. dir,east.dir as 2000 


reserve AQ. of. north. unit ,AQ. of. south. unit,AQ. of. east. unit, 
AW. of. north. unit ,AW. of. south. unit,AW. of. east. unit as 2000 
END ‘'INITIALIZE 


EVENT NORTH. UNIT. GENERATOR saving the event notice 
fej - Ni or js < Ni or je < Ni 
schedule a new. north. unit now 
reschedule this north. unit. generator in 
gamma. f(rl*alphal, betal,sl) minutes 
always 
END ‘' NORTH. UNIT. GENERATOR 


EVENT SOUTH. UNIT. GENERATOR saving the event notice 
if jn < Nl or js < N1 or je < Nil 
schedule a new. south. unit now 
reschedule this south. unit. generator in 
gamma. f(rl*alphal, betal,sl) minutes 
always 
END '' SOUTH. UNIT. GENERATOR 


EVENT EAST. UNIT. GENERATOR saving the event notice 
Jf jn < Nieee js <#Nl*oreje< No 
schedule a new. east. unit now 
reschedule this east.unit. generator in 
gamma. f(r3*alphal,betal,sl) minutes 
always 
END ''EAST. UNIT. GENERATOR 
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EVENT NEW. NORTH. UNIT 
define U as an rea! variable 
create an unit 
let in = Site 
current. time = time. v*hours. v*minutes. v 
arrival. time. of. north. unit(in) = cusneneeeane 
U = uniform. £(0,1,s3) 


if north. run <> 0 or n.north. queue > 0 
go to ‘case3' 


else 
if seed 
next. north. route = 1 
go to ‘'casel' 
else 
next.north. route = 2 
go to 'case2' 
‘casel' 
if Case. runt <> i 
north. run = l 


next. north. route = 0 
jn = jn+il 
schedule an north.service giving unit 

in gamma. f(alpha2,beta2,s2) minutes 
departure. time. of.north. unit( jn) = current. time 
AQ. of. north. unit( jn) = AQ. north. unit 


else 
‘case3' 
file this unit in the north. queue 
north. queue. length2(in) = n. north. queue 
always 
go to ‘end' 
‘case2' 


if south. run = 0 and east.run <> ] 
let north. run = 2 
jn = jnvt+il 
next.north. route = 0 
schedule an north.service giving unit 
in gamma. f(alpha2,beta2,s2) minutes 
departure. time. of. north. unit( jn) = current. time 
AQ. of. north. unit( jn) = AQ.north. unit 
else 
go to ‘case3' 
always 


end’ 

print 1 line with current. time,U, 
north. run, Souch, .un seas nuns. 
n.north. queue,n. south. queue,n. east. queue, 
next. north. route,next. south. route, 


next. east. route thus 
NNU dete dete _ wv sive keke seat seve ese otek Sek tetk sesk 


END " "NEW. NORTH. UNIT 
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EVENT NEW. SOUTH. UNIT 


define U as real variable 

create an unit 

is =is +1 

current. time = time. v*hours. v*minutes. v 
arrival. time. of. south. unit(is) = current. time 
let U = uniform. £(0,1,s3) 


if n.south. queue <> 0 or south. run <> 0 
g0 to ‘case3' 
else 
Hee Oe 1 
next. south. route = 
go to ‘casel' 
else 
next. south. route 
go to ‘case2' 


| 
— 


II 
ho 


casel' 

if north. run <> 2 and east. run = 0 
south. run = 1 
next. south. route = 0 
Hoo 7S] 


schedule an south.service giving unit 

in gamma. f(alpha2,beta2,s2) minutes 
departure. time. of. south. unit( js) = current. time 
AQ. of. south. unit( js) = AQ. south. unit 


else 
case3' 
file this unit in the south. queue 
south. queue. length2(is) = n. south. queue 
always 
go to ‘end' 
case2' 


Peanorel, mim <> 2 
let south. run = 2 
js = js tl 
next. south. route = 0 
schedule an south. service giving unit 
in gamma. f(alpha2,beta2,s2) minutes 
departure. time. of. south. unit( js) = current. time 
AQ. of. south. unit( js) = AQ. south. unit 
else 
go to ‘case3' 
always 


end' 

print 1 line with current. time,U, 
north. run,south. run,east. run, 
n. north. queue,n. south. queue,n. east. queue, 
next. north. route,next. south. route, 


next. east. route thus 
NSU deve deve , vere Seis Gis SBSs eae) Coe wae wes Perle seve 


12 


END ''NEW. SOUTH. UNIT 


EVENT NEW. EAST. UNIT 
define U as real variable 
create an unit 
ie = ie + 1 
current. time = time. v*hours. v*minutes. v 
arrival.time. of.east.unit(ie) = current. time 
U = uniform. £(0,1,s3) 


if n.east. queue > 0 or east. run <> 0 
go to ‘case3' 


else 

if U < Pl 
next. east. route = ]l 
go to ‘casel' 

else 


next.east. route = 2 
go to ‘case2' 
'casel' 
if south. run <> 1 and north. run = 0 
east.run = 1" 
next. east. route = 0 
je = je til 
schedule an east.service giving unit 
in gamma. f(alpha2,beta2,s2) minutes 
departure. time. of. east. unit( je) = current. time 
AQ. of. east. unit( je) = AQ. east. unit 


else 
'case3' 
file this unit in the east. queue 
east. queue. length2(ie) = n. east. queue 
always 
go to ‘end’ 
'case2' 


if soueh rune. 
east.run = 2 
je = jet+l 
next. east. route = 0 
schedule an east.service giving unit 
in gamma. f(alpha2,beta2,s2) minutes 
departure. time. of. east. unit( je) = current. time 
AQ.of.east.unit(je) = AQ. east. unit 
else 
go to ‘case3' 
always 
‘end’ 
print 1 line with current. time,U, 
north. run,south. run,east. run, 


n.north. queue,n. south. queue,n. east. queue, 


next. north. route,next. south. route, 
next. east. route thus 


Be 


NEU 


END 


deve ded , weve Yeve ede deve deve deve deve Fede deve deve 
' "NEW. EAST. UNIT 


EVENT NORTH. SERVICE given unit 

define unit as a pointer variable 

define U as real variable 

destroy this unit 

kn = kn + 1 

current. time = time. v*hours. v*minutes. v 

end. time. of. north. service(kn) = current. time 
te uniform. f(0,1,s3) 


if n.north. queue > 0 and U < Pl 


next. north. route = 1 


always 


if n.north. queue > 0 and U >= Pl 


next. north. route = 2 


always 


if north. run = 2 


else 


if south. run <> 1 and next. east. route 


' ° ! 
go to service2 


north. dir(kn) = 1 
north. run = 0 


[I 
bh 


east. run = l 
je = je +1 
departure. time. of. east. unit(je) = current. time 
waiting. time. of.east.unit(je) = 

current.time - arrival. time. of. east. unit(je) 
east. queue. lengthi(je) = n. east. queue 
remove the first unit from the east. queue 
AQ. of. east.unit( je) = AQ. east. unit 
next. east. route = 0 
schedule an east.service giving unit 

in gamma. f(alpha2,beta2,s2) minutes 


always 


‘next’ 


if east.run <> 1 and next. north. route = 1 


jn = jnt+il 
north. run = 1 
departure. time. of. north. unit( jn) 
waiting. time. of. north. unit( jn) 

current.time - arrival. time. of. north. unit( jn) 
north. queue. lengthl( jn) = n.north. queue 
remove the first unit from the north. queue 
AQ. of. north. unit( jn) = AQ. north. unit 
next. north. route = 0 
schedule an north.service giving unit 

in gamma. f(alpha2,beta2,s2) minutes 


current. time 


always 
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if east. run <> 1 and next. north. route = 2 and south. run = 0 
jn = jntil 
north. run = 2 
departure. time. of. north. unit( jn) = 
waiting. time. of. north. unit( jn) = 
current.time - arrival. time. of. north. unit( jn) 
north. queue. lengthi( jn) = n.north. queue 
remove the first unit from the north. queue 
AQ. of. north. unit( jn) = AQ. north. unit 
next. north. route = 0 
schedule an north. service giving unit 
in gamma. f(alpha2,beta2,s2) minutes 


current. time 


always 
go to ‘print’ 


"service2' 

north. dir(kn) = 2 
north. run = 0 

if next. east. route = 1 


je = jetil 
east. run = l 
departure. time. of. east. unit( je) = current. time 
waiting. time. of. east. unit(je) = 

current. time - arrival. time. of. east. unit( je) 
east. queue. lengthl( je) = n. east. queue 
remove the first unit from the east. queue 
AQ. of. east. unit( je) = AQ. east. unit 
next. east. route = 0 
schedule an east. service giving unit 

in gamma. f(alpha2,beta2,s2) minutes 
always 


if east. run = 0 and next. south. route = 1 
js =js til 
south. run = 1 
departure. time. of. south. unit( js) = current. time 
waiting. time. of. south. unit( js) = current.time - 
arrival. time. of. south. unit( js) 
south. queue. lengthl1(js) = n. south. queue 
remove the first unit from the south. queue 
AQ. of. south. unit( js) = AQ. south. unit 
next. south. route = 0 
schedule an south.service giving unit 
in gamma. f(alpha2,beta2,s2) minutes 
always 
if next. south. route = 2 
js =js tl 
south. run = 2 
departure. time. of. south. unit( js) = current. time 
waiting. time. of. south. unit(js) = 
current. time - arrival. time. of. south. unit(js) 
south. queue. lengthi( js) = n. south. queue 
remove the first unit from the south. queue 
AQ. of. south. unit( js) = AQ. south. unit 
next. south. route = 0 
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schedule an south. service giving unit 
in gamma. f(alpha2,beta2,s2) minutes 
always 
go to ‘next’ 


sorint ' 
print 1 line with current. time,U, 
north. run,south. run,east. run, 
n. north. queue,n. south. queue,n. east. queue, 
next. north. route,next. south. route, 
next. east. route thus 
NS teak deve _ eek Wee IIe | «(IKI ae dbbs cok iG Se eikhs ck 


END ''NORTH. SERVICE 


EVENT SOUTH. SERVICE given unit 
define U as real variable 
define unit as a pointer variable 
ks = ks + 1 
destroy this unit 
current. time = time. v*hours. v*minutes. v 
end. time. of. south. service(ks) = current. time 
U = uniform. f£(0,1,s3) 


if n.south. queue > 0 and U < Pl 
next. south. route = 1 
always 


if n.south. queue > 0 and U >= Pl 
next. south. route = 2 
always 


if south. run = 2 
go to ‘service2’ 
else 
south. dir(ks) 
south. run 


1 
0 


if next. east. route = 1 and north. run = 0 
je = je + il 
east. run = 1 
departure. time. of. east. unit( je) = current. time 
waiting. time. of. east. unit( je) = 
current. time - arrival. time. of. east. unit( je) 
east. queue. lengthl( je) = n. east. queue 
remove the first unit from the east. queue 
AQ. of. east.unit( je) = AQ.east.unit 
next. east. route = 0 
schedule an east.service giving unit 
in gamma. f(alpha2,beta2,s2) minutes 
always 


if next. east. route = 2 
je = je +l 
east. run = 2 
departure. tire. of. east. unit( je) = current. time 


i2 


waiting. time. of. east. unit( je) = current.time - 
arrival. time. of. east. unit( je) 
east. queue. lengthl( je) = n. east. queue 
remove the first unit from the east. queue 
AQ. of. east.unit(je) = AQ. east. unit 
next. east. route = 0 
schedule an east.service giving unit 
in gamma. f(alpha2,beta2,s2) minutes 
always 
"next! 
if next.north. route = 2 and east.run <> 1 
jn = jn+il 
north. run = 2 
departure. time. of. north. unit( jn) = current. time 
waiting. time. of. north.unit( jn) = 
current. time - arrival. time. of. north. unit( jn) 
north. queue. lengthl( jn) = n. north. queue 
remove the first unit from the north. queue 
AQ. of. north. unit( jn) = AQ. north. unit 
next. north. route = 0 
schedule an north.service giving unit 
in gamma. f(alpha2,beta2,s2) minutes 
always 


if east. run = O and north. run <> 2 and next. south. route = 1 
js =jst+l1 . 
south. run = l 
departure. time. of. south. unit( js) = current. time 
waiting. time. of. south. unit(js) = 
current. time - arrival. time. of. south. unit( js) 
south. queue. lengthl( js) = n. south. queue 
remove the first unit from the south. queue 
AQ. of. south. unit( js) = AQ. south. unit 
next. south. route = 0 
schedule an south.service giving unit 
in gamma. f(alpha2,beta2,s2) minutes 
always 


II 
ho 


if north. run <> 2 and next. south. route 
js = jst+l 
south. run = 2 
departure. time. of. south. unit( js) = current. time 
waiting. time. of. south. unit( js) = 
current.time - arrival. time. of. south. unit( js) 
south. queue. lengthl( js) = n. south. queue 
remove the first unit from the south. queue 
AQ. of. south. unit( js) = AQ. south. unit 
next. south. route = 0 
schedule an south. service giving unit 
in gamma. f(alpha2,beta2,s2) minutes 
always 
zo to | print’ 


"service2' 
south. dir(ks) = 2 
south. run = 0 
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go to ‘next’ 


[peant 
eermt 1 line with current. tame,U, 
north. run,south. run,east. run, 
n.north. queue,n. south. queue,n. east. queue, 
next.north. route,next. south. route, 
next.east. route thus 
SS deve dete , ied ise i ga cisel se rhacis Neve Weak aeee wistis Cisgisl  Fiseis 


END ''SOUTH. SERVICE 


EVENT EAST. SERVICE given unit 

define U as real variable 

define unit as a pointer variable 

destroy this unit 

ke = ke + l 

current.time = time. v*hours. v*minutes. v 
end. time. of. east. service(ke) = current. time 
U = uniform. f(0,1,s3) 


if the east.queue is not empty and U < Pl 
next. east. route = l 
always 


if the east. queue is not empty and U >= Pl 
next. east. route = 2 
always 


if east. run = 2 
go to ‘service2' 


else 
east. dir(ke) = 1 
east. run = 0 
if next.north. route = 1 


jn = jnertil 
north. run = 1 
departure. time. of. north. unit( jn) = current. time 
waiting. time. of.north. unit( jn) = 

CULTeNE tamieu- Aalraval amesot norch.unit( jn) 
north. queue. lengthi( jn) = n. north. queue 
remove the first unit from the north. queue 
AQ. of. north. unit( jn) = AQ. north. unit 
next. north. route = 0 
schedule an north.service giving unit 

in gamma. f(alpha2,beta2,s2) minutes 
always 


if next.north. route = 2 and south. run <> 2 
june jn 1 
north. run = 2 
departure. time. of. north. unit( jn) = current. time 
waiting. time.of.north.unit(jn) = 
current.time - arrival. time.of.north. unit( jn) 
north. queue. lengthl( jn) = n.north. queue 


8] 


remove the first unit from the north. queue 
AQ. of. north. unit( jn) = AQinortheunie 
next. north. route = 0 
schedule an north.service giving unit 
in gamma. f(alpha2,beta2,s2) minutes 
always 


‘next’ 
if next. south. route = 1 and north. run <> 2 
js = js + 1 
south. run = 1 
departure. time. of. south. unit( js) = current. time 
waiting. time. of. south. unit( js) = 
current.time - arrival. time. of. south. unit( js) 
south. queue. lengthl( js) = n. south. queue 
remove the first unit from the south. queue 
AQ. of. south. unit( js) = AQ. south. unit 
next. south. route = 0 
schedule an south.service giving unit 
in gamma. f(alpha2,beta2,s2) minutes 
always 


if north. run = 0 and south. run <> 1 and next. east. route = 1 
je =jert+l 
east. run =1 
departure. time. of. east. unit( je) = current. time 
waiting. time. of. east. unit(je) = 
current. time - arrival. time. of. east. unit(je) 
east. queue. lengthl( je) = n. east. queue 
remove the first unit from the east. queue 
AQ. of. east.unit( je) = AQ. east. unit 
next. east. route = 0 
schedule an east.service giving unit 
in gamma. f(alpha2,beta2,s2) minutes 
always 


if south. run <> 1 and next. east. route = 2 
je = je til 
east.run = 2 
departure. time. of. east. unit( je) = current. time 
waiting. time. of. east.unit( je) = 
current.time - arrival. time. of. east. unit(je) 
east. queue. lengthli( je) = n. east. queue 
remove the first unit from the east. queue 
AQ. of.east.unit( je) = AQ. east. unit 
next. east. route = 0 
schedule an east.service giving unit 
in gamma. f(alpha2,beta2,s2) minutes 
always 
26. tO) apaimes 


"service2' 
east. dir(ke) = 2 
east. run = OQ 


go to ‘next’ 
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otent | 


meant 1 line with current. time,U, 


ES 


HOorunarun, SOUCh  Tulmeds te. FUN, 

n. north. queue,n. south. queue,n. east. queue, 
next. north. route,next. south. route, 

next. east. route thus 


, weds sevc ovevs eve weve skeve deve seve sede seve 


END ''EAST. SERVICE 


ROUTINE to SINGLE. RUN. REPORT 
define nn,ns,ne as integer variables 


print 3 line thus 





NO 


ror nn 


do 








Single. run. report — 
WT oT ET Q2 DIR 


SS SS a Ee ee ee ee ee ee Ee ee eee 








print 1 line with nn, 
aerival. time. of. north. unat( nn)» 

north. queue. length2(nn), 

waiting. time. of. north. unit(nn), 
departure. time. of. north. unit(nn), 

end. time. of. north. service(nn), 

poe th. gueue. ee cue nE) north. dir(nn) thus 


slententa 
rer 


loop 


for ns 
do 


s'e s"s . sf e sf A s'- ate : sles ale we ve seve ve ve se se 


print 1 line with ns, 

arrival. time. of. south. unit(ns), 
south. queue. length2(ns), 

waiting. time. of. south. unit(ns), 
departure. time. of. south. unit(ns), 
end. time. of. south. service(ns), 


sesese 


loop 


for ne 
do 


sooth. geuse Fe ita sQucl: daetus) puus 


print 1 line with ne, 
arrival. time. of. east. unit(ne), 

east. queue. length2(ne), 
waiting. time. of. east. unit(ne), 
departure. time. of. east. unit(ne), 

end. time. of. east. service(ne), 

east. queue. lengthl(ne),east.dir(ne) thus 


veveve 
loop 


weds dove deve deve vede | deve seve seve 


END ''SINGLE. RUN. REPORT 
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ROUTINE 
define 
define 


define 


to) SUMSRERORE 

ii as integer variables 

SUM. Of MOT tae win, 

sum. of. south. wt, 

sum. of.east.wt as 1l-dimensional real array 
sum. north. wt, 

sum. south. wt, 

sum.east.wt as real variables 


reserve sum. of.north.wt, 


for ii 
do 


sum. of. south. wt, 
sum. of. east. wt as 2000 


= 1 to Nl 


sum. north. wt = sum.north.wt + 


waiting. time. of. north. unit(ii) 


sum.of.north.wt(ii) = sum. north. wt 
AW. of. north. unit( ii) = sum. of. north. wt(ii)/ii 


Sum. south. wt = sum. south. wt + 


waiting. time. of. south. unit(ii) 


sum. of. south. wt(ii) = sum. south. wt 
AW. of. south. unit(ii) = sum. of. south. wt(ii)/ii 


Sum. east.wt = sum. east.wt + 


waiting. time. of. east. unit(ii) 


sum. of. east.wt(ii) = sum. east. wt - 
AW. of. east. unit(ii) = sum. of. east.wt(ii)/ii 


‘print running average of waiting time and queue length. 


if ii=1 
Prints saline cus 


==SSS= Sum. report routine ===S>===== 


AQn AWn AQs AWs AQe AWe 








Sm rcs cr mc crc rm mm mm ee ee 


PLANE 


1 line with ii, 
AQ. of. north. unit( ia) JAW. of tect. unit Cie 
AQ. of. south. unit( ii) ,AW. of. south. unit(ii), 
AQ. of. east. unit(ii), 
AW. of. east. unit(ii) thus 


desks teve deve tek deve wx vex doe dese deve aede Tere deve 


loop 


release sum. of.north. wt,sum. of. south. wt,sum. of. east. wt 


Stop | 


END SUM. REPORT 
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APPENDIX C. A PROGRAM LISTING OF MODEL 2. 


EVENT NEW. NORTH. UNIT 
define U as an real variable 
create an unit 
let in = in+1 
current. time = time. v*hours. v*minutes. v 
arrival. time. of. north. unit(in) = current. time 
De— uniform. £(0,1,s3) 


if north. run <> 0 or n.north. queue > 0 
go to 'case3' 
else 
shat Ul Ss deral 
next. north. route 
go to ‘casel' 
else 
next. north. route = 2 
go to ‘case2' 


I 
— 


‘casel' 
if east. run <> 1 
north. run = 1 
next. north. route = 0 
Jn = gn 1 
schedule an north.service giving unit 
in gamma. f(alpha2,beta2,s2) minutes 
departure. time. of. north. unit( jn) = current. time 
AQ. of. north. unit( jn) = AQ. north. unit 


else 
‘case3' 
file this unit in the north. queue 
north. queue. length2(in) = n. north. queue 
always 
go to ‘end’ 
‘case2' 


if south. run <> 2 
let north. run = 2 
jn = jn +1 
next.north. route = 0 
schedule an north.service giving unit 
in gamma. f(alpha2,beta2,s2) minutes 
departure. time. of. north. unit( jn) = current. time 
AQ. of. north. unit( jn) = AQ. north. unit 
else 
go to ‘case3' 
always 


end' 
print 1 line with current. time,JU, 
north. run,south. run,east. run, 
n. north. queue,n. south. queue,n. east. queue, 
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next. north. route ,next. south. route, 


next. east. route thus 
NNU wie ede _ ert seve ese | ses P\stis  mie)e iis seve sete | ses 


END ' "NEW. NORTH. UNIT 


EVENT NEW. SOUTH. UNIT 
define U as real variable 
create an unit 
is = is + 1] 
current. time = time. v*hours. v*minutes. v 
arrival. time. of. south. unit(is) = current. time 
let U = uniform. f(0,1,s3) 


if n. south. queue <> 0 or south. run <> 0 
go to ‘case3' 
else 
ot gee URE en) 
next. south. route 
go to ‘casel' 
else 
next. south. route 
go to ‘case2' 


II 
— 


II 
No 


‘Cased. 

if east. run <> 2 
south. run = 1 
next. south. route = 0 
js = js til 


schedule an south.service giving unit 
in gamma. f(alpha2,beta2,s2) minutes 
departure. time. of. south. unit( js) = current. time 
AQ. of. south. unit( js) = AQ. south. unit 
else 


‘case3' 
file this unit in the south. queue 
south. queue. length2(is) = n. south. queue 
always 
go to ‘end’ 
"case2' 
if north. rune 22 
let south. run = 2 
js = js +l 
next. south. route = 0 
schedule an south. service giving unit 
in gamma. f(alpha2,beta2,s2) minutes 
departure. time. of. south. unit( js) = current. time 
AQ. of. south. unit(js) = AQ. south. unit 
else 
Gomtom cases. 
always 
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‘end’ 

print 1 line with current. time,U, 
north. run,south. run,east. run, 
n. north. queue,n. south. queue,n. east. queue, 
next. north. route,next. south. route, 


next. east. route thus 
NSU *x*, *#%* _ ek Soe ie eee Seve oes wor 2 cisely Gk’ aeay 


END ''NEW. SOUTH. UNIT 


EVENT NEW. EAST. UNIT 
define U as real variable 
create an unit 
ie = ie +1 
current. time = time. v*hours. v*minutes. v 
arrival. time. of. east. unit(ie) = current. time 
Oe— uniform. f(0,1,s3) 


if n.east. queue > 0 or east. run <> Q 
go to ‘case3' 
else 
5g ge Oa ame op 
next.east. route = 1 
comEa cases 
else 
next. east. route = 2 
Somber case? 
'casel' 
meenorti.run <> 1 
east.run = 1 
next. east. route = 0 
je = je + 1 
schedule an east.service giving unit 
in gamma. f(alpha2,beta2,s2) minutes 
departure. time. of. east. unit( je) = current. time 
AQ. of. east. unit(je) = AQ. east. unit 
else 
"ease3' 
file this unit in the east. queue 
east. queue. length2(ie) = n. east. queue 
always 
go to ‘end’ 


"case2' 
Peesouth. run <> a 
east.run = 2 
je = je til 
next. east. route = 0 
schedule an east.service giving unit 
in gamma. f(alpha2,beta2,s2) minutes 
departure. time. of. east.unit( je) = current. time 
AQ. of. east. unit( je) = AQ. east. unit 
else 
go to 'case3' 
always 
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t t 


end 
print 1 line with currneme, cme, U, 
north. run,south. run,east. run, 
n. north. queue,n. south. queue,n. east. queue, 
next. north. route,next. south. route, 


next. east. route thus 
NEU 9 2%, vex _ ink Yerke deve deve ‘eve desk deve deve tek ded 


END ''NEW. EAST. UNIT 


EVENT NORTH. SERVICE given unit 

define unit as a pointer variable 

define U as real variable 

destroy this unit 

kn = kn + l 

current. time = time. v*hours. v*minutes. v 

end. time. of. north. service(kn) = current. time 
U = uniform. £(0,1,s3) 


if n.north. queue > 0 and U < Pl 
next. north. route = 1 
always 


if n.north. queue > 0 and U >= Pl 
next. north. route = 2 
always 


if north. run = 2 
go Lor service? 
else 
north. dir(kn) = 1 
north. run = 0 


if next. east. route = 1 
east. run = 1 
je = jetil 
departure. time. of. east. unit( je) = current. time 
waiting. time. of.east.unit(je) = 
current.time - arrival. time. of. east. unit(je) 
east. queue. lengthl( je) = n. east. queue 
remove the first unit from the east. queue 
AQ. of. east. unit( je) = AQ. east. unit 
next. east. route = 0 
schedule an east.service giving unit 
in gamma. f(alpha2,beta2,s2) minutes 
always 


"next' 

if east. run <> 1 and next.north. route = 1 
jn = jntil 
north. run = l 
departure. time. of.north. unit( jn) = 
waiting. time.of.north.unit(jn) = 

current. time - arrival. time. of. north. unit( jn) 

north. queue. lengthi( jn) = n. north. queue 


current. time 
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always 


remove the first unit from the north. queue 
AQ. of. north. unit( jn) = AQ.north. unit 
next. north. route = 0 
schedule an north. service giving unit 
in gamma. f(alpha2,beta2,s2) minutes 


if next. north. route = 2 and south. run <> 2 


jn = jnt+il 
north. run = 2 
departure. time. of.north. unit( jn) = current. time 
waiting. time. of.north. unit(jn) = 

current. time - arrival. time. of. north. unit( jn) 
north. queue. lengthl( jn) = n.north. queue 
remove the first unit from the north. queue 
AQ. of. north. unit( jn) = AQ. north. unit 
next. north. route = 0 
schedule an north.service giving unit 

in gamma. f(alpha2,beta2,s2) minutes 


always 

Posto print 
"service2' 

north. dir(kn) = 2 
Horen. run = 0 


if next. south. route = 2 


js = js + il 
south. run = 2 
departure. time. of. south. unit( js) = current. time 
waiting. time. of. south. unit( js) = 

current.time - arrival. time. of. south. unit( js) 
south. queue. lengthl( js) = n. south. queue 
remove the first unit from the south. queue 
AQ. of. south. unit( js) = AQ. south. unit 
next. south. route = 0 
schedule an south. service giving unit 

in gamma. f(alpha2,beta2,s2) minutes 


always 

go to ‘next’ 

‘print’ 

Brint 1 line with current. time,U, 


NS 


north. run,south. run,east. run, 
n. north. queue,n. south. queue,n. east. queue, 
next. north. route,next. south. route, 


next. east. route thus 
Were vee , ire sieve desde obse Feve seske seve gore Sese coke 


END ' 'NORTH. SERVICE 


EVENT SOUTH. SERVICE given unit 
define U as real variable 
define unit as a pointer variable 
ks = ks + 1 
destroy this unit 
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current. time = time. v*hours. v*minutes. v 
end. time. of. south. service(ks) = current. time 
U = uniform. f(0,1,s3) 


if n. south. queue > 0 and U < Pl 
next. south. route = 1 
always 


if n.south. queue > 0 and U >= Pl 
next. south. route = 2 
always 


if south. run = 2 
go to 'service2' 
else 
south. dir(ks) 
south. run 


1 
0 


if next. east. route = 2 

je =jetil 

east. run = 2 

departure. time. of. east. unit( je) = current. time 

waiting. time. of.east.unit( je) = current.time - 

arrival. time. of. east. unit( je) 

east. queue. lengthl( je) = n. east. queue 

remove the first unit from~-the east. queue 

AQ. of. east.unit( je) = AQ. east. unit 

next. east. route = 0 

schedule an east.service giving unit 

in gamma. f(alpha2,beta2,s2) minutes 

always 


‘next’ 
if east. run <> 2 and next. south. route = 1 
{Se>_ iS oceeel 
south. run = ] 
departure. time. of. south. unit( js) = current. time 
waiting. time. of. south. unit(js) = 
current. time - arrival. time. of. south. unit( js) 
south. queue. lengthl( js) = n. south. queue 
remove the first unit from the south. queue 
AQ. of. south. unit( js) = AQ. south. unit 
next. south. route = 0 
schedule an south.service giving unit 
in gamma. f(alpha2,beta2,s2) minutes 
always 


Il 
ho 


if north. run <> 2 and next. south. route 
js =jst+l 
south. run = 2 . 
departure. time. of. south. unit( js) = current. time 
WATtdIe Gimeno L- SOUCH ale) ina ee alo 

current.time - arrival. time. of. south. unit( js) 

south. queue. lengthl( js) = n. south. queue 
remove the first unit from the south. queue 
AQ. of. south. unit( js) = AQ. south. unit 
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next. south. route = 0 
schedule an south.service giving unit 
in gamma. f(alpha2,beta2,s2) minutes 
always 
go to ‘print’ 


"service2' 
south. dir(ks) = 2 
South. run = 0 


if next.north. route = 2 and south.run <> 2 
fn — jn + 1 
north. run = 2 
departure. time. of. north. unit( jn) = current. time 
waiting. time.of.north.unit(jn) = 
current.time - arrival. time. of.north. unit( jn) 
north. queue. lengthl( jn) = n. north. queue 
remove the first unit from the north. queue 
AQ. of. north. vnit( jn) = AQ. north. unit 
next.north. route = 0 
schedule an north.service giving unit 
in gamma. f(alpha2,beta2,s2) minutes 
always 
go to ‘next’ 


‘print’ 
print 1 line with current.time,U, 
north. run,south. run,east. run, 
n.north. queue,n. south. queue,n. east. queue, 
next. north. route,next. south. route, 
next. east. route thus 
SS vou’ | eve c* Seve Sere Gree seve seve Fess seve seve Vive 


END ‘SOUTH. SERVICE 


EVENT EAST. SERVICE given unit 
define U as real variable 

define unit as a pointer variable 
destroy this unit 


ke = ke + 1 
current. time = time. v*hours. v*minutes. v 
end. time. of. east. service(ke) = current. time 


Oe uniform. £(0,1,s3) 


if the east.queue is not empty and U < Pl 
next.east. route = 1 
always 


if the east.queue is not empty and U >= Pl 
next. east. route = 2 
always 


if east. run = 2 

go to 'service2' 
else 

east.dir(ke) = 1 
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if next. north. route 


always 


t 
"next 


il 
© 


east. run 


ll 
-_ 


ce be 
north. run = 1 
departure. time. of. north. unit( jn) = current. time 
waiting. time.of.north.unit( jn) = 

current.time - arrival. time. of. north. unit( jn) 
north. queue. lengthl( jn) = n. north. queue 
remove the first unit from the north. queue 
AQ. of. north. unit( jn) = AQ. north. unit 
next. north. route = 0 
schedule an north.service giving unit 

in gamma. f(alpha2,beta2,s2) minutes 


if north. run <> 1 and next. east. route = l 


always 


je = je +l 
east. run =1 
departure. time. of. east. unit( je) = current. time 
waiting. time. of. east. unit(je) = 

current. time - arrival. time. of. east. unit( je) 
east. queue. lengthl( je) = n. east. queue 
remove the first unit from the east. queue 
AQ. of. east. unit(je) = AQ. east. unit 
next. east. route = 0 
schedule an east.service giving unit 

in gamma. f(alpha2,beta2,s2) minutes 


if south. run <> 1 and next. east. route = 2 


always 
go to 


je = jetil 
east. run = 2 
departure. time. of. east. unit( je) = current. time 
waiting. time. of.east.unit(je) = 

current.time - arrival. time. of. east. unit(je) 
east. queue. lengthl( je) = n. east. queue 
remove the first unit from the east. queue 
AQ.of.east.unit(je) = AQ. east. unit 
next. east. route = 0 
schedule an east.service giving unit 

in gamma. f(alpha2,beta2,s2) minutes 


Vom ante 


'service2' 
east. dir(ke) = 2 
east. run = 0 


if next. south. route = 1 


js =jst+il 
south. run = 1 
departure. time. of. south. unit( js) = current. time 
waiting. time. of. south. unit(js) = 
current. time - arrival. time. of. south. unit( js) 


a2 


south. queue. lengthl( js) = n. south. queue 

remove the first unit from the south. queue 

AQ. of. south. unit( js) = AQ. south. unit 

next. south. route = 0 

schedule an south.service giving unit 

in gamma. f(alpha2,beta2,s2) minutes 

always 
go to ‘next’ 


‘print’ 

print 1 line with current. time,U, 
north. run,south. run,east. run, 
n. north. queue,n. south. queue,n. east. queue, 
next. north. route,next. south. route, 


next. east. route thus 
ES Tee Het ik die akie sede wewe ved aed wei eke aed 


END ''EAST. SERVICE 
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